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SECTION  A 


EVIDENCE  FOR  Rg-TO-P  SCATTERING  IN  TELESEISMIC  P-CODA  AND 

i 

Rg-TO-S  SCATTERING  IN  REGIONAL  Lg  OF  EAST  KAZAKH  EXPLOSIONS 


SUMMARY 


Analysis  of  teleseismic  data  from  both  Degelen  and  Shagan  River  (East  Kazakh, 
U.S.S.R.)  underground  nuclear  explosions  recorded  at  two  seismic  arrays  suggests  that  the 
low-frequency  teleseismic  P-coda  is  dominated  by  near-source  scattering  of  fundamental-mode 
Rayleigh-to-P  (/?g->P).  The  data  set  consisted  of  22  Degelen  (including  8  with  announced 
shot  depths)  and  40  Shagan  River  explosions  recorded  at  the  NORSAR  and  EKA  arrays, 
respectively.  Good  agreement  between  the  observed  and  expected  variation  of  the  amplitude 
ratio  P/P-coda  (averaged  over  the  frequency  range  of  0.5-2.0  Hz)  with  shot  depth  indicates 
that  scattering  of  the  explosion-generated  Rg,  apparently  due  to  heterogeneities  in  the  source 
region,  makes  significant  contributions  to  teleseismic  P-coda.  Similarly,  regional  data  from 
Shagan  River  explosions  recorded  at  the  CDSN  station,  WMQ,  indicates  that  the  low- 
frequency  Lg  is  probably  due  to  the  scattering  of  explosion-generated  Rg  into  S.  The  scatter¬ 
ing  mechanisms  provide  possible  explanations  for  the  stability  of  P-coda  and  Lg  as  estimators 
of  explosive  yields  because  the  scattering  is  expected  to  occur  over  a  large  region  of  the  shal¬ 
low  crust.  An  understanding  of  the  generation  of  teleseismic  P-coda  and  low-frequency  Lg 
from  underground  nuclear  explosions  not  only  provides  useful  source  information  including 
shot  depth  but  is  also  essential  for  improving  yield  estimates  and  source  discrimination. 
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SECTION  B 


GENERATION  OF  LOW-FREQUENCY  Lg  FROM  NTS  EXPLOSIONS 

SUMMARY 

Analysis  of  regional  data  from  Yucca  Flat  and  Pahute  Mesa  regions  of  NTS  suggests  that 
near-source  scattering  makes  a  significant  contribution  to  the  large-amplitude  low-frequency 
Lg  observed  at  regional  distances.  A  possible  mechanism  for  generation  of  low-frequency  Lg 
appears  to  be  the  scattering  of  explosion-generated  Rg  into  S.  Variations  of  low-frequency 
amplitude  ratios  Pn/Lg  and  Pg/Lg  with  shot  depth  are  not  always  consistent  with  the  expected 
depth-dependence  of  Rg  generated  by  a  point  source  within  a  homogeneous  half  space.  The 
lack  of  good  agreement  in  certain  cases  may,  however,  be  due  to  the  complexity  of  sub¬ 
surface  geology  at  the  NTS  where  large  lateral  and  vertical  variations  in  physical  properties  of 
rock  are  common.  For  closely  spaced  shots  at  both  Yucca  Flat  and  Pahute  Mesa,  a  com¬ 
parison  of  the  low-frequency  spectra  of  their  regional  phases  at  common  recording  stations 
shows  variation  with  shot  depth  that  agrees  well  with  the  theoretically  expected  dependence 
on  shot  depth  and  therefore  supports  the  Rg-to-S  mechanism.  The  generation  of  Lg  appears 
to  be  considerably  influenced  by  azimuthal  direction,  perhaps  because  the  scattering  function 
Rg-to-S  varies  significantly  with  direction.  The  observed  strong  dependence  on  azimuth 
appears  to  support  near-source  scattering  rather  than  spall  as  the  main  contributor  to  the  low- 
frequency  Lg.  Observed  differences  in  the  Lg  spectra  of  normal  and  overburied  shots  may 
also  be  explained  in  terms  of  Rg-to-S  scattering.  The  known  stability  of  Lg  as  an  estimator  of 
explosive  yields  may  be  due  to  scattering  occurring  over  a  relatively  large  volume  of  the 
uppermost  crust. 


IV 


SECTION  C 

AN  ARRAY  STUDY  OF  THE  EFFECTS  OF  A  KNOWN  LOCAL  SCATTERER 

ON  REGIONAL  PHASES 


SUMMARY 


Our  f-k  analyses  of  NORESS  array  recordings  of  telcseismic  events  first  indicated  a 
prominent  secondary  source  about  25-30  km  southwest  of  the  array,  in  the  region  of  Lake 
Mjosa  with  large  topographical  relief.  The  use  of  residual  seismograms  derived  by  subtract¬ 
ing  the  beamed  record  from  each  array  channel  provided  an  effective  method  for  isolating  the 
secondary  source  from  the  effects  of  the  primary  source.  In  this  study,  similar  f-k  techniques 
are  used  to  investigate  the  effects  of  this  known  scatterer  on  the  regional  phases  Pn,  Pg,  and 
Lg  from  several  events  at  regional  distances.  The  beamed  record  for  each  regional  phase  is 
obtained  by  incorporating  steering  delays  based  on  the  known  azimuthal  direction  and  phase 
velocity  from  the  f-k  analysis  of  raw  data.  Both  explosion  and  earthquake  sources  along  vari¬ 
ous  azimuthal  directions  to  NORESS  provide  evidence  for  scattered  arrivals  Pn-to-Rg,  Pg-to- 
Rg,  and  Lg-to-Rg,  originating  from  the  same  region  of  Lake  Mjosa.  These  scattered  arrivals 
are  observed  with  remarkable  clarity  and  closeness  to  the  expected  arrival  times  and  are  rich 
in  low  frequencies,  as  expected  for  Rg.  Their  presence  in  various  regional  phases  indicates 
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SECTION  A 


EVIDENCE  FOR  Rg-TO-P  SCATTERING  IN  TELESEISMlC  P-CODA  AND 
Rg-TO-S  SCATTERING  IN  REGIONAL  Lg  OF  EAST  KAZAKH  EXPLOSIONS 

INTRODUCTION 

Theoretical  (e.g.  Garvin,  1956)  and  observational  studies  (e.g.  Rulev,  1965)  indicate 
that,  in  the  vicinity  of  a  shallow  explosion,  the  largest  amount  of  ground  surface  motion  is 
generally  due  to  Rayleigh  waves.  An  analysis  of  teleseismic  P  arrivals  from  underground 
nuclear  explosions  at  Novaya  Zemlya,  U.S.S.R.,  known  for  its  rugged  topography,  led 
Greenfield  (1971)  to  suggest  that  a  significant  portion  of  arrivals  in  P-coda  (taken  to  be 
between  10  and  40  sec  after  first  P)  may  be  due  to  near-surface  scattering  of  explosion¬ 
generated  Rayleigh  waves  into  teleseismic  P.  Gupta  et  al.  (1985a,  1991b)  proposed  that 
near-source  scattering  of  fundamental-mode  Rayleigh  waves  (Rg)  into  teleseismic  P  (referred 
to  as  Rg->P  hereafter)  may  be  responsible  for  the  large  arrivals  immediately  following  P  and 
pP  on  NORSAR  records  of  Nevada  Test  Site  (NTS)  shots.  Theoretical  studies  of  near-source 
scattering  by  explosion  sources  in  the  Yucca  Flat  (NTS)  region  by  McLaughlin  et  al.  (1987) 
and  Stead  and  Helmbergcr  (1988)  also  demonstrated  that  locally  scattered  Rayleigh  waves 
may  be  responsible  for  the  large  arrivals  immediately  following  P  and  pP.  F-k  analyses  of 
teleseismic  recordings  of  both  U.S.  and  U.S.S.R.  explosions  provided  evidence  in  favor  of 
both  near-receiver  P-»Rg  and  near-source  Rg-*P  scattering  (Gupta  et  al. ,  1990a;  Gupta 
et  al. ,  1990b).  Observations  of  the  spectra  of  Rg  signals  from  both  explosions  and  very 
shallow-focus  earthquakes  indicate  the  dominant  energy  in  Rg  to  be  confined  to  frequencies 
less  than  about  2  Hz  and  the  most  important  parameter  for  influencing  its  generation  to  be  the 
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source-depth  (e:g..- Kafka,  1990). 


Several  investigators  (e.g.  Frasier,  1972;  Bakun  and  Johnson,  1973;  McLaughlin  et  ai, 
199.0)' have  suggested  that  spall  makes  a  significant  contribution  to  the  low-frequency  com¬ 
ponent  of  teleseismic  P  arrivals.  It  is  important  to  understand  the  generation  of  teleseismic 
Prcoda,  .especially  because  of  its  usefulness  as  an  effective  measure  of  the  seismic  yield  of 
nuclear  ^explosions;  (Gupta  et.al.,  1985b).  Moreover,  anunderstanding  of  the  various  arrivals 
;in  .teleseismic  !P  from  nuclear  explosions  may  provide  useful  source  information  such  as  shot 
depth,  medium  velocity,  ,  and  .geological  environment.  In  this  study,  we  analyze  the  spectra  of 
teleseismic  ,‘P  and  :;P-coda  from  underground  nuclear  explosions  at  'Degelen  Mountain  and 
:Shagan  iRiyer  test  sites  in  Eastern  ‘Kazakh,  U.S.S.R.,  recorded  at  the  NORSAR  and  EKA 
arrays,  respectively.  Least  squares  inversion  of  array  data  allows  :the  results  to  be  relatively 
free  from  recording  site  effects.  Use  is  made  of  the  tecently  released  information  on  Soviet 
nuqlear  explosions  by  Bocharov  et  al.  (1989).  A  comparison  of  the  observed  dependence  of 
the  amplitude  ratio  P/P-coda  on  shot  depth  with  the  theory  of  explosion-generated  Rg  shows 
good  agreement,  suggesting  that  the  low-frequency  content  of  teleseismic  P-coda  consists 
mainly  of  Rg-»P  arrivals.  Similarly,  analysis  of  regional  data  from  14  Shagap  River  explo¬ 
sions  recorded  at  the  CDSN  station,  WMQ,  indicates  that  the  low-frequency  component  of  the 
observed  Lg  is  probably  due  to  the  scattering  of  explosion-generated  Rg  into  S.  It  is  expected 
that  Rg-*P  scattering  should  be  accompanied  by  even  larger  Rg->S  scattering  (e.g.  Stead  and 
Hefinberger,  1988). 
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DEPENDENCE  OF  EXPLOSION-GENERATED  Rg  AND  P-CODA  ON  SHOT  DEPTH 


For  an  explosion  source  with  source  function  S(f)  within  a  homogeneous  half  space,  the 
fundamental-mode  Rayleigh  wave  amplitude  Rg(f)  is  given  by  the  expression  (Hudson  and 
Douglas,  1975): 


10 


-4.4  —  f 


where  <j)(co)  is  the  Fourier  transform  of  the  source  RDP  <j>(t),  h  is  the  shot  depth,  to  =  2ttf,  a  is 
the  compressional-wave  velocity  of  the  shot  medium,  p  is  the  density,  r  is  the  epicentral  dis¬ 
tance,  A  is  a  numerical  constant,  and  a  Poisson’s  ratio  of  0.25  is  assumed.  The  far-field  P- 
wave  displacement,  in  the  frequency  domain,  may  be  written  as  (Hudson  and  Douglas,  1975): 


P(f)  = - —  to  <J)(0))  (2) 

4  7t  p  cr  R 

where  R  is  the  teleseismic  source-receiver  distance.  If  the  source-receiver  path  and  the 
medium  velocity  do  not  change,  equations  (1)  and  (2)  yield 


(f)  =  k  C(f)  f  °-5  10  ”4’4  a  f 


where  k  is  a  constant  and  C(f)  denotes  the  scattering  function  Rg->P.  Assuming  that  the 
low-frequency  P-coda  in  teleseismic  data  is  mainly  due  to  the  scattering  of  Rg,  we  replace 
Rg— »P  in  equation  (3)  by  P-coda  and  obtain 

log  — (0  =  4.4  ±  f  -  log  k  -  log  C(f)  -  0.5  log  f  (4) 

P-coda  a 

This  shows  that  if  we  consider  explosions  for  which  the  Rg-»P  scattering  function  may  be 
assumed  to  be  the  same  (such  as  closely  spaced  explosions)  and  the  frequency  is  kept  fixed 
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(i.c.  for  a  narrow  frequency  band),  a  plot  of  log  P/P-coda  versus  shot  depth,  h  should  have  a 
slope  of  about  4.4  f/a. 
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P  AND  P-CODA  FROM  DEGELEN  EXPLOSIONS  RECORDED  AT  NORSAR  ARRAY 

We  first  examined  short-period,  vertical  component  telescismic  P  arrivals  from  22 
Degelen  explosions  (Table  1)  recorded  at  22  center  element  sensors  of  the  NORSAR  array, 
providing  over  400  observations.  In  Table  1,  data  for  the  first  8  shots  is  from  Bocharov  et  al. 
(1989)  and  those  for  the  remaining  explosions  is  from  P.  D.  Marshall  (written  communication, 
1987).  The  epicentral  distance  is  about  38°  along  the  azimuth  of  about  313°  or  N  47°  W. 
This  testing  area  is  dominated  by  the  Degelen  mountain,  with  local  relief  of  a  few  hundred 
meters  and  locally  steep  slopes.  Figure  1,  based  on  topographic  data  in  Bonham  et  al. 
(1980),  shows  the  main  features  along  with  the  eight  known  locations  of  shots  from  Bocharov 
et  al.  (1989).  The  emplacement  of  nuclear  devices  was  made  possible  by  tunneling  into  the 
mountain. 

As  compared  to  results  from  a  single  receiver,  use  of  multichannel  data  has  the  advan¬ 
tage  that  receiver  effects,  which  can  drastically  influence  the  observed  signals,  can  be  elim¬ 
inated  or  at  least  suppressed.  Excluding  a  few  channels  with  spikes  or  clipped  data,  spectra 
and  spectral  ratios  P/P-coda  were  obtained  for  each  record  by  selecting  a  window  length  of 
12.8  sec  (beginning  4  sec  before  the  onset  of  P)  as  the  P-window  and  the  following  25.6  sec 
signal  as  the  P-coda,  and  applying  Parzen  taper  to  each  (see  Figure  1  of  Gupta  and  Blandford, 
1987  for  examples).  A  correction  for  noise  was  made  by  selecting  a  sample  of  noise  before 
the  onset  of  P,  obtaining  its  spectra,  and  subtracting  its  power  from  that  of  the  observed  sig¬ 
nal.  The  average  amplitude  ratio  of  P/P-coda,  over  the  frequency  range  of  0.5-2.0  Hz,  was 
obtained  by  using  only  those  data  points  for  which  the  signal/noise  power  ratio  was  at  least  2 
for  each  of  P  and  P-coda.  All  available  amplitude  ratios  P/P-coda  were  used  as  input  to  a 
least  squares  inversion  that  separated  the  source  and  receiver  terms. 
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Figure  1.  Topographic  map  of  the  Degelen  mountain  region  with  locations  of  8  explosions 
numbered  as  in  Table  L  The  contour  interval  is  200  m. 


TABLE  1 


22  DEGELEN  MOUNTAIN  SHOTS  RECORDED  AT  NORSAR  ARRAY 


No. 

DATE 

ORIGIN  TIME 

LAT  (N) 

LON  (E) 

mb 

DEPTH 

(M) 

YIELD 

(KT) 

1 

22  Mar  1971 

04:33:00.3 

49.79847 

78.10897 

5.767 

283.0 

20-150 

2 

25  May  1971 

04:03:00.4 

49.80164 

78.13883 

5.048 

132.0 

<20 

3 

29  Nov  1971 

06:02:59.9 

49.74342 

78.07850 

5.462 

203.0 

<20 

4 

30  Dec  1971 

06:21:00.2 

49.76003 

78.03714 

5.838 

249.0 

20-150 

5 

10  Mar  1972 

04:56:59.8 

49.74531 

78.11969 

5.453 

171.0 

<20 

6 

28  Mar  1972 

04:22:00.1 

49.73306 

78.07569 

5.177 

124.0 

6 

7 

07  Jun  1972 

01:28:00.0 

49.82675 

78.11547 

5.422 

208.0 

20-150 

8 

16  Aug  1972 

03:16:59.8 

49.76547 

78.05883 

5.105 

139.0 

8 

9 

16  Feb  1973 

05:02:57.6 

49.822 

78.158 

5.477 

- 

- 

10 

26  Oct  1973 

04:26:57.7 

49.759 

78.164 

5.228 

- 

- 

11 

16  May  1974 

03:02:57.7 

49.752 

78.093 

5.242 

- 

- 

12 

25  Jun  1974 

03:56:57.7 

49.840 

78.166 

4.419 

- 

- 

13 

10  Jul  1974 

02:56:57.6 

49.783 

78.130 

5.168 

- 

- 

14 

13  Sep  1974 

03:02:57.6 

49.778 

78.081 

5.144 

- 

- 

15 

16  Dec  1974 

06:40:58.0 

49.867 

78.087 

4.887 

- 

- 

16 

20  Feb  1975 

05:32:57.6 

49.789 

78.062 

5.765 

- 

17 

11  Mar  1975 

05:42:57.6 

49.747 

78.146 

5.424 

- 

- 

18 

07  Aug  1975 

03:56:57.7 

49.812 

78.161 

5.187 

- 

- 

19 

15  Jan  1976 

04:46:57.6 

49.824 

78.201 

5.182 

- 

- 

20 

21  Apr  1976 

04:57:57.8 

49.776 

78.146 

4.940 

- 

- 

21 

19  May  1976 

02:56:57.8 

49.796 

78.058 

4.722 

- 

- 

22 

23  Jul  1976 

02:32:57.8 

49.779 

78.085 

4.961 

- 

- 

A  plot  of  the  source-term  amplitude  ratio  P/P-coda  versus  shot  depth  for  8  Degelen 
explosions  for  which  shot  depths  are  known  (Bocharov  et  al. ,  1989)  is  shown  in  Figure  2a 
whereas  a  plot  of  the  source  term  versus  mb  for  all  22  shots  is  shown  in  Figure  2b.  Assuming 
the  P-wave  velocity  in  the  Degelen  testing  region  to  be  about  4  km/sec  (Bonham  et  al. , 
1980),  4.4  f/a  =  1.4,  if  f  is  taken  to  be  the  average  of  0.5  and  2.0,  or  1.25  Hz.  Good  agree¬ 
ment  between  the  mean  slope  expected  from  simple  theory  and  observation  (Figure  2a)  sug¬ 
gests  that  the  low-frequency  P-coda  is  mainly  due  to  the  scattering  of  explosion-generated  Rg. 
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Figure  2.  Source  term  amplitude  ratio  P/P-coda  (log  units),  averaged  over  the  frequency 
range  of  0.5-2.0  Hz  and  derived  from  NORSAR  array  data,  plotted  versus  (a)  shot  depth  for  8 
Degelen  explosions  (numbered  as  in  Table  1),  and  (b)  versus  mb  for  22  Degelen  explosions. 
The  least  squares  linear  regression  (dashed  line),  mean  slope  (with  associated  standard  devia¬ 
tion),  and  one  standard  deviation  of  residuals  (SD)  are  indicated  on  each  plot. 
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Results  in  Figure  2b,  in  which  the  8  explosions  used  in  Figure  2a  have  been  identified,  sug¬ 
gest  generally  larger  P/P-coda  for  larger  explosions  and  a  dependence  on  shot  depth  that 
appears  to  be  similar  to  that  in  Figure  2a. 

On  the  basis  of  the  known  locations  of  8  Degelen  shots  (Bocharov  et  al.,  1989),  Figure 
2a  includes  results  from  three  pairs  of  explosions  (nos.  1  and  2,  3  and  6,  4  and  8)  that  are 
within  about  2  km  of  each  other.  Since  differences  in  propagation  paths  are  known  to  exert 
significant  influence  on  various  seismic  arrivals,  results  from  closely  spaced  seismic  sources 
recorded  at  a  common  station  should  be  considered  more  reliable  than  others.  It  is  important 
to  note  that  the  slope  from  each  of  the  three  pairs  of  closely  spaced  shots  in  Figure  2a  is 
nearly  the  same  as  that  derived  from  the  average  of  all  (8)  shots  and  is  in  good  agreement 
with  the  theoretically  expected  value. 
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P  AND  P-CODA  FROM  SHAGAN  RIVER  EXPLOSIONS  RECORDED  AT  ERA  ARRAY 


We  examined  the  short-period,  vertical-component  data  from  40  Shagan  River  explosions 
(28  from  SW  and  12  from  NE  region)  recorded  at  the  Eskdalemuir  (ERA)  array.  These  shots 
are  listed  in  Table  2  in  which  data  for  the  first  4  shots  is  from  Bocharov  et  al.  (1989)  and 
those  for  the  remaining  explosions  is  from  P.  D.  Marshall  (written  communication,  1987). 
The  array  has  20  elements  arranged  in  ah  "L"  configuration  with  element  spacing  of  0.9  km 
and  maximum  spacing  of  9.8  km  (Bache  et  al. ,  1985).  The  epicentral  distance  of  Shagan 
explosions  to  the  array  is  about  47°  along  the  azimuth  of  about  310°  or  N50°W.  The  Shagan 
River  test  site  is  believed  to  be  underlain  by  about  100  m  of  unconsolidated  sediments  with 
the  thickness  decreasing  gradually  towards  southwest  (Bonham  etal.,  1980).  A  plot  of  the 
source-term  amplitude  ratio  P/P-coda  versus  mb  (Figure  3a)  suggests  generally  larger  P/P-coda 
for  larger  explosions  and  the  mean  slope  is  close  to  that  obtained  iti  Figure  2b.  For  a  fixed 
value  of  mb,  the  NE  Shagan  shots  in  Figure  3a  have  generally  smaller  values  of  P/P-coda  than 
the  SW  shots.  This  will  be  consistent  with  generally  larger  P-coda  and  greater  complexity  for 
the  NE  shots  than  for  the  SW  shots  observed  by  Bache  etal.  (1984)  and  Marshall  et  al. 
(1985). 
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TABLE  2 


40  SHAGAN  RIVER  EXPLOSIONS  RECORDED  AT  EKA  ARRAY 


No. 

DATE 

ORIGIN  TIME 

LAT  (N) 

LON  (E) 

mi> 

DEPTH 

(M) 

YIELD 

(KT) 

1 

15  Jan  1965 

06:00:00.8 

49.93500 

79.00936 

5.87 

178.0 

100-150 

2 

19  Jun  1968 

05:05:59.8 

49.98025 

78.98550 

5.28 

316.0 

<20 

3 

30  Nov  1969 

03:32:59.7 

49.92428 

78.95575 

6.02 

472.0 

125 

4 

30  Jun  1971 

03:56:59.8 

49.94600 

78.98047 

4.94 

217.0 

<20 

5 

31  May  1974 

03:27:00.4 

49.950 

78.852 

5.81 

- 

- 

6 

16  Oct  1974 

06:33:00.5 

49.979 

78.898 

5.41 

- 

- 

7 

27  Dec  1974 

05:46:59.8 

49.943 

79.011 

5.50 

- 

- 

8 

27  Apr  1975 

05:37:00.2 

49.949 

78.926 

5.51 

- 

- 

9 

30  Jun  1975 

03:27:00.5 

50.004 

78.957 

4.52 

- 

- 

10 

29  Oct  1975 

04:47:00.3 

49.946 

78.878 

5.61 

- 

- 

11 

25  Dec  1975 

05:17:00.1 

50.044 

78.814 

5.69 

- 

- 

12 

04  Jul  1976 

02:57:00.4 

49.909 

78.911 

5.85 

- 

- 

13 

28  Aug  1976 

02:57:00.4 

49.969 

78.930 

5.74 

- 

- 

14 

07  Dec  1976 

04:57:00.3 

49.922 

78.846 

5.80 

- 

- 

15 

29  May  1977 

02:57:00.5 

49.937 

78.770 

5.75 

- 

- 

16 

29  Jun  1977 

03:07:00.7 

50.006 

78.869 

5.20 

- 

- 

17 

05  Sep  1977 

03:03:00.3 

50.035 

78.921 

5.73 

- 

- 

18 

30  Nov  1977 

04:07:00.3 

49.958 

78.885 

5.89 

- 

- 

19 

1 1  Jun  1978 

02:57:00.5 

49.898 

78.797 

5.83 

- 

- 

20 

05  Jul  1978 

02:47:00.4 

49.887 

78.871 

5.77 

- 

- 

21 

15  Sep  1978 

02:37:00.3 

49.916 

78.879 

5.89 

- 

- 

22 

04  Nov  1978 

05:06:00.2 

50.034 

78.943 

5.56 

- 

- 

23 

23  Jun  1979 

02:57:00.4 

49.903 

78.855 

6.16 

- 

- 

24 

07  Jul  1979 

03:47:00.2 

50.026 

78.991 

5.84 

- 

25 

04  Aug  1979 

03:57:00.0 

49.894 

78.904 

6.13 

- 

- 

26 

18  Aug  1979 

02:52:00.0 

49.943 

78.938 

6.13 

- 

- 

27 

28  Oct  1979 

03:16:59.9 

49.973 

78.997 

5.98 

- 

- 

28 

02  Dec  1979 

04:37:00.4 

49.891 

78.796 

5.99 

- 

- 

29 

23  Dec  1979 

04:57:00.4 

49.916 

78.755 

6.13 

- 

- 

30 

29  Jun  1980 

02:33:00.6 

49.939 

78.815 

5.69 

- 

- 

31 

12  Oct  1980 

03:34:17.0 

49.961 

79.028 

5.88 

- 

- 

32 

14  Dec  1980 

03:47:09.3 

49.899 

78.938 

5.93 

- 

- 

33 

27  Dec  1980 

04:09:11.0 

50.057 

78.981 

5.87 

- 

- 

34 

29  Mar  1981 

04:03:52.9 

50.007 

78.982 

5.49 

- 

- 

35 

22  Apr  1981 

01:17:14.3 

49.885 

78.810 

5.94 

- 

- 

36 

13  Sep  1981 

02:17:21.2 

49.910 

78.915 

6.06 

- 

- 

37 

18  Oct  1981 

03:57:05.6 

49.923 

78.859 

6.00 

- 

- 

38 

29  Nov  1981 

03:35:11.5 

49.887 

78.860 

5.62 

- 

« 

39 

31  Aug  1982 

01:31:03.6 

49.924 

78.761 

5.20 

- 

. 

40 

26  Dec  1982 

03:35:17.1 

50.071 

78.988 

5.58 

- 
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Figure  3.  Source  term  amplitude  ratio  P/P-coda  (log  units),  averaged  over  the  frequency 
range  of  0.5-2.0  Hz,  for  40  Shagan  River  explosions  recorded  at  the  EKA  array  plotted  versus 
(a)  mb  and  (b)  estimated  shot  depth.  The  least  squares  linear  regression  (dashed  line),  mean 
slope  (with  associated  standard  deviation),  and  one  standard  deviation  of  residuals  (SD)  are 
indicated  on  each  plot.  Six  pairs  of  closely  spaced  explosions  are  identified  by  the  letters  A 
through  F  and  the  SW  and  NE  explosion  populations  are  denoted  by  X  and  O,  respectively. 
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A  relationship  between  the  average  ratio  P/P-coda  and  shot  depth  was  obtained  by  con¬ 
verting  the  known  mb  into  shot  depth  by  using  empirical  relationships  based  on  the  informa¬ 
tion  on  shot  depths  and  yields  recently  released  by  Bocharov  et  al.  (1989).  These  relation¬ 
ships,  derived  from  the  Shagan  River  and  Konystan  (Murzhik)  explosions  are  (R.  S.  Jih,  writ¬ 
ten  communications,  1990): 

mb  =  0.690  log  Y  +  4.605  (5) 

where  Y  is  the  yield  in  kt,  and 

log  h  =  0.241  log  Y  +  2.174  (6) 

where  h  is  shot  depth  in  meters.  These  two  equations  would  indicate  an  explosion  of  mb  = 
6.0  to  have  a  yield  of  about  105  kt  and  a  depth  of  458  m.  Using  equations  (5)  and  (6),  one 
obtains 

log  h  =  0.349  mb  +  0.567  (7) 

which  can  be  used  to  estimate  shot  depth  for  an  explosion  with  known  mb.  A  plot  of  the 
amplitude  ratio  P/P-coda  versus  estimated  shot  depth  in  km  (Figure  3b)  indicates  a  mean  slope 
of  0.943  with  one  standard  deviation  of  residuals  of  only  0.095.  Assuming  the  P-wave  velo¬ 
city  for  the  uppermost  0.5  km  of  the  crust  in  the  Shagan  test  site  region  to  be  5  km/sec,  4.4 
f/a  =  1.1,  if  f  is  taken  to  be  the  average  of  0.5  and  2.0,  or  1.25  Hz.  Good  agreement  between 
the  mean  slope  expected  from  simple  theory  and  observation  (Figure  3b)  suggests  that  the 
low-frequency  P-coda  is  mainly  due  to  the  scattering  of  explosion-generated  Rg. 

Figures  3a  and  3b  include  results  from  six  (including  two  from  NE  Shagan)  pairs  of 
explosions  within  less  than  2  km  of  each  other  and  with  their  mb  differing  by  at  least  0.4 
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m.u.;  these  are  identified  by  the  letters  A  through  F.  Note  that  the  slope  from  each  of  the  six 
pairs  of  shots  in  Figure  3b  is  nearly  the  same  as  derived  from  the  average  of  all  (40)  shots 
and  is  in  good  agreement  with  the  theoretically  expected  value. 

Results  similar  to  those  in  Figure  3  but  for  the  frequency  range  of  3.0-5.0  Hz  are  shown 
in  Figure  4  which  shows  good  separation  between  explosions  from  the  SW  and  NE  regions  of 

5 

the  Shagan  River  test  site.  The  lack  of  any  systematic  variation  with  shot  depth  is  consistent 
with  the  suggestion  of  explosion-generated  Rg  contributing  to  P-coda  because  Rg  is  important 
only  for  frequencies  less  than  about  2  Hz. 
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Figure  4  Similar  to  Figure  3  but  for  the  frequency  range  of  3.0-5. 0  Hz.  There  is  no 
significant  dependence  of  P/P-coda  on  mb  or  estimated  shot  depth  but  the  SW  and  NE  explo¬ 
sion  populations  appear  well  separated. 
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Pn  AND  Lg  FROM  SHAGAN  RIVER  EXPLOSIONS  RECORDED  AT  WMQ 


We  examined  the  spectra  of  the  regional  phases  Pn  and  Lg  from  Shagan  River  nuclear 
explosions  in  order  to  understand  their  dependence  on  source  depth  and  crustal  structure 
(Gupta  et  al. ,  199.1a).  The  data  set  included  14  shots  from  the  southwest  region  of  the  test 
site  recorded  on  the  broadband  instrument  at  the  Chinese  network  station  Urumchi  (WMQ), 
located  about  950  km  southeast  of  the  test  site.  The  explosions  are  listed  in  Table  3,  which 
includes  data  from  Ringdal  and  Marshall  (1989)  for  12  explosions  prior  to  the  year  1989. 
The  first  arrival  Pn  was  distinct  on  all  records,  and  the  beginning  of  Lg  was  assumed  to  be  at 
a  group  velocity  of  3.5  km/sec,  in  agreement  with  the  observed  data.  Spectra  of  Pn  and  Lg 
were  obtained  by  applying  a  10%  cosine  taper  to  time  windows  of  12.8  and  51.2  sec  for  Pn 
aqd  Lg,  respectively.  Average  Pn/Lg  amplitude  ratios,  with  correction  for  noise  and  no 
smoothing,  were  computed  for  desired  ranges  of  frequency. 


TABLE  3 


14  SHAGAN  RIVER  (SOUTHWEST)  EXPLOSIONS  RECORDED  AT  WMQ 


No. 

DATE 

ORIGIN  TIME 

LAT  (N) 

LON  (E) 

mb 

1 

12  Mar  1987 

01:57:17.2 

49.939 

78.823 

5.31 

2 

03  Apr  1987 

01:17:08.0 

49.928 

78.829 

6.12 

3 

20  Jun  1987 

00:53:04.8 

49.913 

78.735 

6.03 

4 

02  Aug  1987 

00:58:06.8 

49.880 

78.917 

5.83 

5 

15  Nov  1987 

03:31:06.7 

49.871 

78.791 

5.98 

6 

13  Dec  1987 

03:21:04.8 

49.989 

78.844 

6.06 

7 

27  Dec  1987 

03:05:04.7 

49.864 

78.758 

6.00 

8 

13  Feb  1988 

03:05:05.9 

49.954 

78.910 

5.97 

9 

03  Apr  1988 

01:33:05.8 

49.917 

78.945 

5.99 

10 

04  May  1988 

00:57:06.8 

49.928 

78.769 

6.09 

11 

14  Sep  1988 

04:00:00.0 

49.870 

78.820 

6.03 

12 

17  Dec  1988 

04:18:06.8 

49.818 

78.910 

5.80 

13 

12  Feb  1989 

04:15:06.8 

49.93 

78.74 

5.9 

14 

08  Jul  1989 

03:46:57.6 

49.87 

78.82 

5.6 

16 


In  order  to  explore  the  possibility  that  the  low-frequency  Lg  from  the  underground 
nuclear  explosions  at  Shagan  River  test  site  is  mainly  due  to  the  scattering  of  explosion¬ 
generated  Rg  into  S,  we  again  compare  common-station  results  from  explosions  with 
significantly  different  shot  depths.  Assuming  the  medium  velocity  to  be  the  same  for  all 
explosions,  equations  (1)  and  (2)  give 

(f)  =  k  c,(o  f w  10  ~*A  “  r  (8) 

where  Cj(f)  denotes  the  scattering  function  Rg— >S.  If  the  low-frequency  Lg  in  regional  data 
is  mainly  due  to  the  scattering  of  Rg— >S,  we  can  replace  Rg— >S  and  P  in  equation  (8)  by  Lg 
and  Pn,  respectively  and  obtain 

log  |i(0  =  4.4-£f- log  k  -  log  C,(0  -  0.5  log  f  (9) 

Thus  if  we  consider  explosions  for  which  the  Rg->S  scattering  function  may  be  assumed  to 
be  the  same  (such  as  closely  spaced  explosions)  and  the  frequency  is  kept  fixed  (i.e.  for  a  nar¬ 
row  frequency  band),  a  plot  of  log  Pn/Lg  versus  shot  depth,  h  should  have  a  slope  of  about 
4.4  f/a. 

Figure  5a  shows  a  plot  of  average  amplitude  ratio  Pn/Lg  over  the  low-frequency  range  of 
0.5-2.0  Hz  versus  mb  for  14  explosions  from  the  southwest  region  of  the  Shagan  River  test 
site  for  which  the  assumption  of  the  same  Rg-»S  scattering  function  may  be  approximately 
valid.  Shot  depths  for  these  Shagan  River  explosions  may  again  be  estimated  by  using  equa¬ 
tion  (7).  A  plot  of  the  amplitude  ratio  Pn/Lg  versus  estimated  shot  depth  (Figure  5b)  shows  a 
mean  slope  of  about  1.08,  in  excellent  agreement  with  the  slope  of  about  1.1  expected  from 
theory. 
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Figure  5.  Average  amplitude  ratio  Pn/Lg  over  the  frequency  range  of  0.5-2.0  Hz  versus  (a) 
mb  and  (b)  estimated  shot  depth  for  14  explosions  from  the  southwest  Shagan  test  site.  Data 
from  two  closely  spaced  shots  (nos.  1  and  2)  are  identified.  Each  plot  indicates  the  least 
squares  linear  regression  (dashed  line),  mean  slope  (with  associated  standard  deviation),  and 
one  standard  deviation  of  residuals  (SD). 
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Figures  5a  and  5b  include  results  from  two  explosions  (nos.  1  and  2)  that  arc  within  a 
few  km  of  each  other.  Regional  phases  are  extremely  sensitive  to  differences  in  propagation 
paths  (Blandford,  1981).  In  a  study  of  regional  data  from  earthquakes  and  explosions  in  the 
Western  United  States,  Lynnes  et  al.  (1990)  found  that  for  frequencies  up  to  4  Hz,  the 
source-station  path  is  more  important  than  the  source  type  in  determining  the  Pg/Lg  ratios. 
Therefore,  results  from  seismic  sources  with  nearly  common  paths  should  be  given  more 
weight  than  others.  Note  that  the  slope  from  the  two  closely  spaced  shots  in  Figure  5b  is 
nearly  the  same  as  that  derived  from  the  average  of  all  (14)  shots  and  is  in  good  agreement 
with  the  theoretically  expected  value.  It  seems  therefore  that  the  low-frequency  Lg  is  gen¬ 
erated  by  the  scattering  of  explosion-generated  Rg  into  S. 
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CONCLUSIONS 


Analysis  of  a  large  amount  of  teleseismic  and  regional  data  from  East  Kazakh  under¬ 
ground  nuclear  explosions  suggests  that  the  low-frequency  teleseismic  P-coda  and  Lg  are 
dominated  by  the  near-source  scattering  of  fundamental-mode  Rayleigh  (Rg)  to  P  and  S, 
respectively.  The  amplitude  ratios  P/P-coda  and  Pn/Lg,  averaged  over  the  frequency  range  of 
0.5-2.0  Hz  and  derived  from  teleseismic  and  regional  data,  respectively,  vary  with  shot  depths 
with  a  mean  slope  in  good  agreement  with  that  expected  from  theory.  Scattering  occurring 
over  a  large  region  of  the  uppermost  crust  may  also  contribute  to  the  observed  stability  of  P- 
coda  and  Lg  over  teleseismic  P  and  regional  Pn,  respectively.  Results  of  this  study  not  only 
provide  new  insight  into  the  generation  of  teleseismic  P-coda  and  low-frequency  Lg  but  also 
improved  methods  for  estimating  shot  depths  of  underground  nuclear  explosions.  It  will  be 
useful  to  apply  these  techniques  to  more  teleseismic  and  regional  data  and  estimate  shot 
depths  of  U.S.S.R.  explosions.  Similar  analysis  should  also  be  carried  out  for  NTS  explosions 
for  which  the  near-source  environmental  conditions  (including  shot  depth)  are  known  but  the 
shallow  crustal  structure  is  considerably  more  complex  than  in  East  Ka/.akh. 
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SECTION  B 


GENERATION  OF  LOW-FREQUENCY  Lg  FROM  NTS  EXPLOSIONS 

INTRODUCTION 

At  regional  distances,  the  S-wave  group  known  as  Lg  is  often  the  largest  amplitude 
arrival  from  both  earthquake  and  explosion  sources.  .Since  Lg  has  proven  useful  for  detection, 
source  discrimination,  and  yield  estimation  of  underground  nuclear  explosions,  it  is  important 
to  understand  the  generation,  propagation,  and  spectral  characteristics  of  Lg.  Perhaps  the 
most  puzzling  aspect  of  Lg  from  explosions  is  their  spectra’s  relative  richness  in  low- 
frequency  content  compared  to  corresponding  earthquake  Lg  spectra,  at  least  for  the  Nevada 
Test  Site  (NTS)  shots  (Murphy  and  Bennett,  1982).  Gupta  and  Blandford  (1983)  suggested 
that  the  low-frequency  Lg  from  explosions  may  be  due  to  preferential  scattering  of  low  fre¬ 
quencies  in  P  to  S  scattering  near  the  source.  Theoretical  studies  of  Lg  generation  (e.g.  Cam- 
pillo  et  al. ,  1984)  indicate  that  explosions  do  not  generate  large  Lg  amplitudes  unless  they  are 
located  close  to  strong  interfaces  such  as  the  free  surface.  There  are  several  other  possible 
mechanisms  that  can  generate  large  S  waves  and  low-frequency  Lg;  these  include  spall  (Tay¬ 
lor  and  Randall,  1989;  McLaughlin  et  al. ,  1990)  and  interaction  of  the  spherical  P-wavefront 
emanating  from  an  explosion  with  the  free-surface,  resulting  in  shear  waves  referred  to  as  S* 
(Gutowski  et  al.,  1984;  Lilwall,  1988).  Pg  is  also  a  prominent  phase  from  both  explosion  and 
earthquake  sources,  especially  in  the  WUS,  and  its  generation  is  also  not  fully  understood. 
Murphy  and  Bennett  (1982)  found  the  Pg  spectra  also  to  be  deficient  in  high-frequency  energy 
as  compared  to  those  from  earthquakes  whereas  the  Pn  spectra  for  the  two  types  of  sources 
were  found  to  be  indistinguishable. 
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As  shown  in  Part  A  of  this  report,  our  analysis  of  regional  data  from  Shagan  River 
explosions  recorded  at  the  CDSN  station,  WMQ,  indicates  that  the  low-frequency  Lg  is  prob¬ 
ably  due  to  the  scattering  of  explosion-generated  Rg  into  S.  In  this  study,  we  examine  the 
possibility  that  a  similar  mechanism  is  responsible  for  generating  the  large  low-frequency  Lg 
observed  from  NTS  explosions  for  which  the  near-source  parameters  (including  shot  depth) 
are  known  but  the  shallow  crustal  structure  is  considerably  more  complex  than  in  East 
Kazakh.  Regional  data  from  66  NTS  explosions  recorded  mostly  at  two  broadband  digital 
stations  at  regional  distances  are  analyzed  for  their  dependence  on  shot  depth  and  other  geo¬ 
logical  and  geophysical  parameters. 


LATERAL  VARIATIONS  IN  MEDIUM  VELOCITY  AND  OTHER  PARAMETERS  AT  NTS 

Extreme  lateral  and  vertical  variations  in  medium  velocity  exist  at  both  Yucca  Flat  and 
Pahute  Mesa  test  regions  of  the  NTS.  At  Yucca  Flat,  the  large  variations  are  due  to  the  pres¬ 
ence  of  low  velocity  alluvium  and  tuffs  within  a  basin  structure  and  the  presence  of  Yucca 
Fault  and  other  geological  features  (Hays  and  Murphy,  1971).  Variations  in  medium  velocity 
arc  even  more  dramatic  at  Pahute  Mesa  because  of  the  presence  of  alternating  layers  of  vari- 

t 

ous  tuffs  and  rhyolites  of  irregular  thickness.  In  a  study  of  geophysical  properties  of  shot 
media  at  NTS,  Ramspott  and  Howard  (1975)  observed  that,  "of  the  various  media,  Yucca  Flat 
and  Pahute  Mesa  above  the  water  table  are  the  most  variable.  In  these  areas,  it  is  possible  to 
find  individual  past  sites  with  extreme  values  of  reported  parameters.  The  other  media  are 
relatively  uniform."  It  is  interesting  to  note  that  nearly  all  explosions  with  large  low- 
frequency  (0.5- 1.0  Hz)  Lg,  used  in  Murphy  and  Bennett’s  (1982)  study,  were  above  the  water 
table  (AWT)  and  shallow  (shot  depth  less  than  500m)  or  within  a  region  where  scattering  is 
expected  to  be  large. 

Table  1  provides  a  chronological  list  of  66  NTS  (including  39  Yucca  Flat  and  27  Pahute 
Mesa)  shots  used  in  this  study.  The  Yucca  Flat  shots  include  13  explosions  detonated  below 
the  water  table  (BWT)  and  26  above  the  water  table  (AWT)  whereas  the  Pahute  Mesa  shots 
include  17  BWT  and  10  AWT  explosions.  The  mb  values  are  maximum-likelihood  estimates 
from  P.  D.  Marshall  (written  communication,  1988)  and  the  medium  velocity  in  the  table  is 
the  overburden  velocity  or  the  average  compressional-wave  velocity  between  the  shot  point 
and  the  surface  as  measured  by  sonic  logs.  Most  explosions  were  recorded  at  both  broadband 
digital  stations  ELK  and  KNB  operated  by  the  Lawrence  Livermore  National  Laboratory.  Fig¬ 
ures  1  and  2  show  the  geographical  locations  of  these  shots  along  with  their  shot  depths.  It  is 
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29 ‘BASEBALL  564 
30  ALIGOTE  320 

32  PINEA‘D  204 

33  ISLAY  294 

34  TREBBIANO  305 

35  CERNADA  213 

36  PAL’IZ'A  472 

37  TILCI  445 

38  ;ROUSANNE  518 

39  AKAVI  494 

40  JORNADA  640 
43  TENAJA  357 

45  BOUSCHET  564 

47  ATRISCO  640 

48  CERRO  229 

49  BORREGO  564 

50  SEYVAL  366 

51  MANTECA  413 

52  COALORA  274 

54  TORQUOISE  533 

55  ARMADA  265 

56  FAHADA  384 

57  DANABLU  320 

58  SABADO  320 

60  1 1 -Cl  IADO  533 

61  ROMANO  515 

62  TORTUGAS  640 

63  MUNDO  567 

64  CAPROCK  600 
66  CORREO  335  ‘ 


Figure  1.  Location  map  of  39  Yucca  Flat  shots  with  13  BWT  and  26  AWT  explosions 
denoted  by  x  and  o,  respectively.  The  numbers  correspond  to  those  in  Table  1.  Shot  names 
and  depths  (in  meters)  are  indicated  on  the  right  of  the  figure.  Note  several  pairs  of  closely 
spaced  explosions  with  considerably  different  shot  depths. 
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1  ALMENDRO  1064  10  POOL  879  28  SERPA  573 

2  TYBO  765  11  FONDUTTA  633  31  HARZER  637 

3  STILTON  732  12  BACKBEACH  672  41  MOLBO  638 

4  MAST  911  13  PANIR  681  42  HOSTA  640 

5  CAMEMBERT  1311  14  FARM  689  44  GIBNE  570 

6  INLET  819  15  PEPATO  681  46  NEBBIOLO  640 

7  MUENSTER  1452  21  COLWICK  633  53  CABRA  543 

8  CHESHIRE  1167  24  KASH  645  59  CHANCELLOR  625 

9' ESTUARY  868  25  TAFI  680  65  KAPPELI  640 

Figure  2.  Location  map  of  27  Pahute  Mesa  shots  with  17  BWT  and  10  AWT  explosions 
denoted  by  x  and  o,  respectively.  The  numbers  correspond  to  those  in  Table  1.  Shot  names 
and  depths  (in  meters)  are  indicated  at  the  bottom  of  the  figure.  Note  several  pairs  of  closely 
spaced  explosions  with  considerably  different  shot  depths. 
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interesting:  to  note  that  several  explosions  are  within  a  few.  hundred  meters  of  each  other  so 
that’  propagation  paths  to  a  common  recording:  station  should  be  almost  identical.  Two  such 
examples  im  which  BWF  and:  AWT  shots  are  closely  spaced.’ are  the  Yucca  Flat  shots  PAL1ZA 
(ho:  36)  and'  ATRISCO  (no.  47).  in  Figure  1  and  the  Pahute  Mesa  shots  MAST  (no.  4).  and 
HOSTA  (no.  42)  in  Figure  2. 


TABLE- 1 

66  NTS  EXPLOSIONS  USED  IN  STUDY 


No: 

DATE 

1 

06  Jun  1973 

2 

14  May  1975 

3 

03  Jun  1975 

4 

19  Jun- 1975 

5 

26  Jun  1975 

6 

20  Nov  1975 

7 

03  Jan  1976 

8 

14  Feb  1976 

9 

09  Mar  1976 

10 

17  Mar  1976 

11 

11  Apr  1978 

12 

11  Apr  1978 

13 

31  Aug  1978 

14 

16  Dec  1978 

15 

11  Jun  1979 

16 

20  Jun  1979 

17 

03  Aug  1979 

18 

08  Aug  1979 

19 

06  Sep  1979 

20 

16  Apr  1980 

21 

26  Apr  1980 

22 

02  May  1980 

23 

22  May  1980 

24 

12  Jun  1980 

25 

25  Jul  1980 

26 

31  Jul 1980 

NAME 


ALMENDRO  (E,K) 
TYBO  (K) 

STILTON.  (E,K) 
MAS!'  (E,K) 
CAMEMBERT  (E) 
INLET  (E,K) 
MUENSTER  (E) 
CHESHIRE  (E,K) 
ESTUARY  (E,K) 
POOL  (E,K) 
FONDUTTA  (E,K) 
BACKBEACH  (E,K) 
PANIR  (K) 

FARM  (K) 

PEPATO  (E) 

CHESS  (E) 

BURZET  (E,K) 
OFFSHORE  (E,K) 
HEARTS  (E,K) 
PYRAMID  (E,K) 
COLWICK  (E,K) 
CANFIELD  (K) 
FLORA  (E,K) 

KASH  (K) 

TAFI  (E,K) 
VERDELLO  (E) 


mb 

VEL 

(KM/S) 

6.166 

3.083 

6.046 

2.275 

5.831 

2.150 

6.060 

3.877 

6.212 

2.900 

5.981 

3.042 

6.288 

2.900 

5.942 

3.210 

5.974 

2.660 

6.061 

2.865 

5.499 

2.299 

5.551 

2.043 

5.667 

2.096 
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(KM) 

(KM) 

1.0638 

0.6860 

0.7650 

0.6300 

0.7315 

0.2720 

0.9113 

0.6660 

1.3106 

0.6680 

0.8 190, 

0.7030 

1.4524 

0.6760 

1.1670 

0.6250 

0.8681 

0:6270 

0.8793 

0.6900 

0.6330 

0.6680 

0.6721 

0.7150 

0.6810 

0.6450 

0.6890 

0.6490 

0.6810 

0.5790 

0.3353 
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0.3353 

0.4730 

0.6450 
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0.6800 

0.6070 

0.3658 
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27 

14  Nov  1980 

DAUPHIN  (E,K) 

4.554 

1.420 

0.3200 

28 

17  Dec  1980 

SERPA  (E,K) 

5.294 

2.000 

0.5730 

29 

15  Jan  1981 

BASEBALL  (E,K) 

5.725 

1.970 

0.5639 

30 

29  May  1981 

ALIGOTE  (E,K) 

4.368 

1.400 

0.3200 

31 

06  Jun  1981 

HARZER  (E,K) 

5.625 

2.056 

0.6370 

32 

16  Jul  1981 

PINEAU  (E,K) 

3.290 

1.125 

0.2042 

33 

27  Aug  1981 

ISLAY  (E) 

4.360 

1.440 

0.2940 

34 

04  Sep  1981 

TREBBIANO  (E,K) 

3.820 

1.465 

0.3048 

35 

24  Sep  1981 

CERNADA  (E,K) 

3.500 

1.060 

0.2130 

36 

01  Oct  1981 

PALIZA  (E,K) 

5.115 

1.497 

0.4724 

37 

11  Nov  1981 

TILCI  (E,K) 

5.035 

1.600 

0.4450 

38 

12  Nov  1981 

ROUSANNE  (E,K) 

5.458 

1.580 

0.5182 

39 

03  Dec  1981 

AKA VI  (E,K) 

4.749 

1.730 

0.4940 

40 

28  Jan  1982 

JORNADA  (E,K) 

5.932 

1.695 

0.6400 

41 

12  Feb  1982 

MOLBO  (E,K) 

5.412 

2.005 

0.6380 

42 

12  Feb  1982 

HOSTA  (E,K) 

5.599 

2.020 

0.6401 

43 

17  Apr  1982 

TENAJA  (E,K) 

4.454 

1.310 

0.3570 

44 

25  Apr  1982 

GIBNE  (E,K) 

5.364 

1.920 

0.5700 

45 

07  May  1982 

BOUSCHET  (E,K) 

5.730 

1.740 

0.5639 

46 

24  Jun  1982 

NEBBIOLO  (E,K) 

5.636 

2.027 

0.6400 

47 

05  Aug  1982 

ATRISCO  (E,K) 

5.738 

1.604 

0.6400 

48 

02  Sep  1982 

CERRO  (E,K) 

3.290 

1.165 

0.2286 

49 

29  Sep  1982 

BORREGO  (E,K) 

4.150 

1.860 

0.5639 

50 

12  Nov  1982 

SEYVAL  (K) 

4.545 

1.320 

0.3660 

51 

10  Dec  1982 

MANTECA  (E,K) 

4.741 

1.610 

0.4130 

52 

11  Feb  1983 

COALORA  (E,K) 

4.150 

1.340 

0.2740 

53 

26  Mar  1983 

CABRA  (E,K) 

5.275 

1.757 

0.5430 

54 

14  Apr  1983 

TORQUOISE  (E,K) 

5.742 

1.868 

0.5330 

55 

22  Apr  1983 

ARMADA  (E,K) 

3.977 

1.313 

0.2650 

56 

26  May  1983 

FAHADA  (E,K) 

4.647 

1.500 

0.3840 

57 

09  Jun  1983 

DANABLU  (E) 

4.519 

1.578 

0.3200 

58 

11  Aug  1983 

SABADO  (E) 

4.391 

1.315 

0.3200 

59 

01  Sep  1983 

CHANCELLOR  (E) 

5.445 

1.970 

0.6250 

60 

22  Sep  1983 

TECHADO  (E,K) 

4.150 

1.610 

0.5334 

61 

16  Dec  1983 

ROMANO  (K) 

5.140 

1.699 

0.5150 

62 

01  Mar  1984 

TORTUGAS  (K) 

5.848 

1.769 

0.6400 

63 

01  May  1984 

MUNDO  (K) 

5.474 

1.670 

0.5670 

64 

31  May  1984 

CAPROCK  (K) 

5.743 

1.585 

0.6000 

65 

25  Jul  1984 

KAPPELI  (E,K) 

5.385 

2.100 

0.6400 

66 

02  Aug  1984 

CORREO  (E,K) 

4.669 

1.305 

0.3350 

Note:  Stations  recording  the  event  are  denoted  by  E  (ELK)  and  K  (KNB). 
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Plots  of  delay  time  (=  2h/a  where  h  is  the  shot  depth  and  a  the  overburden  velocity) 
versus  shot  depth  for  39  Yucca  Flat  and  27  Pahute  Mesa  shots  are  shown  in  Figure  3.  The 
Yucca  Flat  data  (Figure  3a)  shows  an  approximately  linear  trend  suggesting  gradual  increase 
in  velocity  with  depth.  But  the  Pahute  Mesa  data  in  Figure  3b  indicate  extreme  variations  in 
medium  velocity  with  depth.  Large  lateral  variations  in  velocity  are  mainly  responsible  for 
the  Jack  of  an  approximately  linear  relationship.  For  example,  the  explosions  MAST  and 
HOSTA,  within  about  500  m  of  each  other  (Figure  2),  have  overburden  velocities  of  about  3.9 
and  2.0  km/sec  so  that  the  delay  time  for  the  deeper  (911  m)  shot,  MAST,  is  smaller  than  that 
for  the  shallower  (640  m)  shot,  HOSTA.  Most  of  these  drastic  variations  are  perhaps  of  lim¬ 
ited  spatial  extent  and  therefore  not  "seen"  by  seismic  waves  of  low  frequency  or  long 
wavelength.  This  means  that,  if  only  low  frequencies  are  being  used,  a  deeper  shot  should 
almost  always  be  considered  to  have  a  larger  delay  time  than  a  shallower  explosion.  In  other 
words,  delay  times  derived  from  the  linear  regression  relationships  in  Figures  3a  and  3b  for 
shots  with  known  depths,  referred  to  hereafter  as  "estimated  delay  times",  may  be  expected  to 
provide  better  estimates  than  those  from  sonic  logs  in  the  immediate  vicinity  of  the  shot  holes. 
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(a)  39  YF  SHOTS 


SHOT  DEPTH  (KM) 


SHOT  DEPTH  (KM) 


Figure  3.  Delay  time  versus  shot  depth  for  (a)  39  Yucca  Flat  and  (b)  27  Pahute  Mesa  shots 
used  in  the  study.  The  least  squares  linear  regression  line,  mean  slope  (with  associated  stan¬ 
dard  deviation),  one  standard  deviation  of  residuals  (SD),  and  the  correlation  coefficient  (CC) 
are  indicated  on  each  plot. 
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COMPARISON  OF  LOW  FREQUENCY  Pn.Pg,  AND  Lg 


We  first  examined  spectra  of  the  regional  phases  Pn,  Pg,  and  Lg  from  a  fairly  large 
number  of  NTS  explosions  recorded  at  two  common  stations,  ELK  and  KNB,  located  at  epi- 
central  distances  of  about  410  and  290  km,  respectively.  In  order  to  isolate  source  effects,  it 
is  important. to  use  common  source-receiver  paths  because  of  the  extreme  sensitivity  to  propa¬ 
gation  paths  of  the  regional  phases,  especially  Pn  (Langston,  1982)  and  Lg  (Gupta  and 
Blandford,  1983).  Source  characteristics  derived  from  closely  spaced  explosions  should  there¬ 
fore  be  given  more  importance  than  others.  The  window  lengths  used  for  spectral  analysis  of 
Pn,  Pg,  and  Lg  on  the  vertical  component  records  were  6.4,  12.8,  and  25.6  sec,  respectively, 
and  a  10%  cosine  taper  was  used.  The  beginning  of  Pn  was  quite  sharp  for  nearly  all  explo¬ 
sions  and  the  onsets  of  Pg  and  Lg  phases  were  selected  in  the  same  manner  as  in  Murphy  and 
Bennett  (1982). 

Amplitude  ratios  Pn/Lg  and  Pg/Lg,  averaged  over  the  frequency  ranges  of  0.3- 1.0  Hz  and 
0.5-2.0  Hz,  were  computed  for  all  explosions  recorded  at  both  ELK  and  KNB  and  plotted 
versus  shot  depth,  delay  time,  and  estimated  delay  time.  In  order  to  avoid  possible  effects  of 
large  variations  in  scaled  depth  on  regional  phases,  7  overburied  shots  (all  from  Yucca  Flat) 
were  excluded  from  plots.  Moreover,  results  from  these  overburied  shots  appeared  as  outliers 
on  most  plots.  Plots  of  all  available  amplitude  ratios  Pn/Lg  and  Pg/Lg  from  the  remaining  59 
shots  versus  shot  depth,  delay  time,  and  estimated  delay  time,  derived  from  ELK  data,  are 
shown  in  Figures  4  through  7.  The  figures  include  results  of  linear  regression  such  as  the 
mean  slope,  standard  deviation  of  residuals  (SD),  and  the  correlation  coefficient  values  (CC). 
Results  from  Yucca  Flat  shots  alone  (Figures  4  and  6)  and  those  from  a  combination  of  Yucca 
Flat  and  Pahute  Mesa  shots  (Figures  5  and  7)  are  shown.  The  scattering  functions  at  Yucca 
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26  YF  ELK  26  YF  ELK 


Figure  4.  Average  amplitude  ratio  Pn/Lg  (log  units)  over  the  frequency  ranges  of  0.3- 1.0  Hz 
and  0.5-2.0  Hz  for  26  Yucca  Flat  shots  recorded  at  ELK  versus  shot  depth,  delay  time,  and 
estimated  delay  time.  Each  plot  indicates  the  least  squares  linear  regression,  mean  slope  (with 
associated  standard  deviation),  one  standard  deviation  of  residuals  (SD),  and  the  correlation 
coefficient  (CC). 
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49  NTS, ELK  Pn6.4/Lg25.6  0.3-1. 0  HZ 


49  NTS  ELK  Pn6.4/Lg25.6  0.5-2.0  HZ 
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Figure  5.  Similar  to  Figure  4  for  49  NTS  (26  Yucca  Flat  and  23  Pahute  Mesa)  shots. 
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Figure  6.  Similar  to  Figure  4  for  Pg/Lg  derived  from  26  Yucca  Flat  shots. 


37 


49  NTS  ELK;Rg12.8/lg25:6  0.3-1 .0  HZ 


49  NTS  ELK  Pg12:8/Lg25.6  0.5-2.0  HZ 


Figure  7.  Similar  to  Figure  4  for  Pg/Lg  derived  from  49  NTS  (26  Yucca  Flat  and  23  Pahute 
Mesa)  shots. 
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Flat  and  Pahute  Mesa  may  not  be  the  same  but  the  combined  results  are  included  so  that  a 
bigger  range  of  source  parameters  can  be  examined. 

The  correlations  in  several  figures,  such  as  Figures  6f  and  7f,  both  for  Pg/Lg  versus 
estimated  delay  time,  are  good' whereas  those  in  Figures  4c  and  5c  are  poor.  In  general, 
linear  relationships  for  the  frequency  range  of  0.5-2.0  Hz  are  better  than  those  for  0.3-1. 0  Hz 
and  those  for  Pg/Lg  are  better  than  those  for  Pn/Lg.  Possible  reasons  for  the  latter  are  greater 
stability  of  Pg  over  Pn,  use  of  signal  time  window  that  is  twice  as  long  as  for  Pn,  and  much 
better  S/N  for  Pg  than  for  Pn  (especially  for  the  low  frequencies).  Theoretical  considerations 
(Gupta  et  al.y  1991)  would  suggest  the  mean  slopes  of  Pn/Lg  versus  delay  time  to  be  about 
2.2  f  where  f  is  the  frequency,  or  about  1.4  for  the  frequency  range  0.3- 1.0  if  f  is  taken  to  be 
the  average  of  0.3  and  1.0,  or  0.65  Hz.  Similarly,  theory  would  suggest  a  mean  slope  of 
about  2.8  for  the  frequency  range  of  0.5-2.0  Hz.  One  may  assume  that  the  depth  dependence 
of  low-frequency  Pg/Lg  is  similar  to  that  for  Pn/Lg  because  Pg  is  derived  mainly  from  P 
waves  trapped  within  the  crust.  Thus  theory  would  suggest  the  mean  slopes  in  Figures  6c  and 
7c  to  be  about  1.4  and  those  in  Figures  6f  and  7f  to  be  about  2.8.  It  should  be  noted  that  the 
mean  slopes  indicated  in  Figures  4  through  7  are  all  based  on  the  assumption  that  values  on 
the  x-axis  (shot  depth,  delay  time,  or  estimated  delay  time)  are  error-free  and  that  all  error  lies 
in  the  measurement  of  the  average  spectral  ratios.  This  assumption  is  unrealistic  for  several 
reasons  such  as  finite  (rather  than  a  point)  source  size  and  large  lateral  variations.  It  is,  how¬ 
ever,  known  that  a  linear  regression  that  incorporates  an  error  in  values  on  the  x-axis  will 
always  increase  the  mean  slope  (see,  for  example,  Bolt,  1978)  and  thus  bring  the  slope  values 
indicated  on  linear  regressions  such  as  those  in  Figures  6c,  6f,  7c,  and  7f  close  to  the  theoreti¬ 
cally  predicted  values. 
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All  available  values  of  amplitude  ratios  Pn/Lg  and  Pg/Lg,  averaged  over  the  frequency 
ranges  of  0.3- 1.0  Hz  and  0.5-2.0  Hz,  are  plotted  at  their  epicentral  (ground  zero)  locations  in 
Figures  8  and  9  for  26  Yucca  Flat  shots  recorded  at  ELK.  Location  plots  of  shot  depth  and 
delay  times  for  the  same  26  explosions  are  included  in  Figure  10.  These  data  include  several 
pairs  of  closely  spaced  explosions  for  which  the  Rg-to-S  scattering  functions  are  expected  to 
be  similar.  Pairs  of  closely  spaced  shots  with  their  depths  differing  by  at  least  50  m  are 
identified  in  Figures  8,  9,  and  i0.  When  the  shallower  shot  is  AWT  and  the  deeper  explosion 
is  BWT,  they  are  denoted  by  the  letters  A  and  B,  respectively.  Other  pairs  are  marked  by  the 
letters  S  and  D  denoting  the  relatively  shallower  and  deeper  shots,  tespectively.  A  com¬ 
parison  of  the  four  plots  in  Figures  8  and  9  with  the  two  plots  in  Figure  10  show  considerable 
similarity.  Moreover,  for  closely  spaced  shots,  the  deeper  shot  (marked  B  or  D)  almost 
always  has  relatively  larger  Pn/Lg  and  Pg/Lg  for  both  ranges  of  frequency  than  those  for  the 
shallower  shot  (denoted  by  A  or  S).  It  seems  therefore  that,  for  closely  spaced  shots,  the 
observed  values  of  Pn/Lg  and  Pg/Lg  are  in  agreement  with  the  theoretically  expected  depen¬ 
dence  on  shot  depth  and  delay  time.  These  observations  also  suggest  that  the  near-source 
Rg-to-S  scattering  functions  are  perhaps  nearly  the  same  only  for  closely  spaced  shots.  The 
lack  of  a  linear  relationship  for  many  plots  in  Figures  4  through  7  may  therefore  be  due  to  the 
scattering  functions  being  significantly  different  from  one  shot  region  to  another,  presumabl) 
due  to  the  lateral  variations  that  exist  at  the  Yucca  Flat  test  site. 

Plots  of  amplitude  ratios  Pn/Lg  and  Pg/Lg,  averaged  over  the  frequencies  ranges  of  0.3- 
1.0  Hz  and  0.5-2.0  Hz  and  obtained  from  23  Pahute  Mesa  shots  recorded  at  ELK,  are  plotted 
at  their  epicentral  locations  in  Figures  11  and  12.  Figure  13  shows  plots  of  shot  depth  and 
delay  time  values  for  the  same  23  Pahute  Mesa  shots.  Unlike  the  remarkable  similarity 
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Figure  8.  Average  amplitude  ratio  Pn/Lg  (in  log  units  and  relative  to  the  indicated  mean 
value)  derived  from  ELK  records  of  26  Yucca  Flat  shots  plotted  at  their  known  epicentral 
locations  for  the  frequency  range  of  (a)  0.3- 1.0  Hz  and  (b)  0.5-2.0  Hz.  Three  pairs  of  closely 
spaced  A  (AWT)and  B  (BWT)  shots  and  two  pairs  of  S  (shallower)  and  D  (deeper)  shots  are 
identified.  For  nearly  all  pairs,  the  deeper  shot  (B  or  D)  has  relatively  larger  Pn/Lg  than  the 


shallower  (A  or  S)  explosion. 
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37.00 


Figure  9.  Similar  to  Figure  8  for  Pg/Lg.  For  nearly  all  pairs,  the  deeper  shot  (B  or  D)  has 
larger  Pg/Lg  than  the  shallower  (A  or  S)  explosion. 
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Figure  10.  Similar  to  Figure  8  for  26  Yucca  Flat  shots  for  (a)  shot  depth  and  (b)  delay  time. 
Note  similarity  of  the  two  figures. 
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Figure  11.  Similar  to  Figure  8  for  Pn/Lg  for  23  Pahute  Mesa  shots.  Six  pairs  of  closely 
spaced  A  and  B  shots  and  one  pair  of  S  and  D  shots  are  identified.  For  nearly  all  pairs,  the 
deeper  shot  (B  or  D)  has  larger  Pn/Lg  than  the  shallower  (A  or  S)  explosion. 
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Figure  13.  Similar  to  Figure  8  for  23  Pahute  Mesa  shots  for  (a)  shot  depth  and  (b)  delay 
time.  Note  the  relative  lack  of  similarity  between  the  two  figures,  especially  for  two  pairs  of 
closely  spaced  explosions  (MAST  and  HOSTA,  INLET  and  BACKBEACH). 


between  plots  of  shot  depth  and  delay  time  for  Yucca  Flat  explosions  (Figures  10a  and  l()b), 
the  Pahute  Mesa  plots  (Figures  13a  and  13b)  show  much  less  similarity.  In  fact,  for  two  pairs 
of  closely  spaced  AWT  and  BWT  explosions  (MAST  and  HOSTA,  INLET  and  BACK- 
BEACH),  the  shot  depth  and  delay  time  have  opposite  trends.  As  mentioned  earlier, 
extremely  large  lateral  variations  in  medium  velocity  are  responsible  for  this.  Figures  11,  12, 
and  13  include  data  from  7  pairs  of  closely  spaced  explosions  such  that  6  pairs  include  AWT 
and  BWT  shots  (denoted  by  the  letters  A  and  B,  respectively)  and  the  remaining  pair  (denoted 
by  S  and  D)  consists  of  two  BWT  shots.  A  comparison  of  the  four  plots  in  Figures  1 1  and 
12  with  the  shot  depth  plot  in  Figure  13a  show  a  great  deal  of  similarity,  especially  for  the 
relative  differences  among  closely  spaced  shots.  As  observed  for  Yucca  Flat  shots,  the  deeper 
shots  (marked  B  and  D)  almost  always  have  relatively  larger  values  of  both  Pn/Lg  and  Pg/Lg 
for  both  ranges  of  frequency  than  those  for  the  shallower  shots  (denoted  by  A  and  S).  Once 
again,  it  seems  that,  for  closely  spaced  shots,  the  observed  values  of  Pn/Lg  and  Pg/Lg  are  in 
agreement  with  the  theoretically  expected  dependence  on  shot  depth  and  the  lack  of  a  linear 
relationship  for  Pahute  Mesa  data  points  in  Figures  5  and  7  are  mainly  due  to  the  strong 
lateral  variations  in  the  geology  of  the  Pahute  Mesa  region. 

Plots  of  all  available  amplitude  ratios  Pn/Lg  and  Pg/Lg  (averaged  over  the  frequency 
ranges  of  0.3- 1.0  Hz  and  0.5-2.0  Hz)  versus  shot  depth,  delay  time,  and  estimated  delay  time, 
derived  from  51  NTS  shots  (including  28  from  Yucca  Flat)  recorded  at  KNB,  are  shown  in 
Figures  14  through  17.  The  correlations  are  poor  in  some  cases  and  good  in  some  others. 
The  results  are  therefore  similar  to  those  for  the  recording  station  ELK  (Figures  4  through  7). 
All  available  values  of  Pn/Lg  and  Pg/Lg  derived  from  KNB  are  plotted  at  their  epicentral 
locations  in  Figures  18,  19,  21,  and  22  whereas  the  corresponding  values  of  shot  depth  and 
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Figure  14.  Similar  to  Figure  4  for  28  Yucca  Flat  shots  recorded  at  KNB. 
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Figure  15.  Similar  to  Figure  4  derived  from  51  NTS  (28  Yucca  Flat  and  23  Pahute  Mesa) 
shots  recorded  at  KNB. 
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Figure  16.  Similar  to  Figure  4  for  Pg/Lg  derived  from  28  Yucca  Flat  shots  recorded  at  KNB. 
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Figure  17.  Similar  to  Figure  4  for  Pg/Lg  derived  from  48  NTS  (28  Yucca  Flat  and  20  Pahute 
Mesa)  shots  recorded  at  KNB. 
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Figure  21.  Similar  to  Figure  8  for  23  Pahute  Mesa  shots  recorded  at  KNB.  Four  pairs  of 
closely  spaced  A  and  B  shots  and  two  pairs  of  S  and  D  shots  are  identified.  For  nearly  all 
pairs,  the  deeper  shot  (B  or  D)  has  relatively  larger  Pn/Lg  than  the  shallower  (A  or  S)  explo¬ 
sion. 
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Figure  22.  Similar  to  Figure  8  for  Pg/Lg  derived  from  20  Pahute  Mesa  shots  recorded  at 
KNB.  For  nearly  all  pairs,  the  deeper  shot  (B)  has  relatively  larger  Pg/Lg  than  the  shallower 
(A)  explosion. 
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Figure  23.  Similar  to  Figure  8  for  23  Pahute  Mesa  shots  recorded  at  KNB  for  (a)  shot  depth 
and  (b)  delay  time.  Note  the  relative  lack  of  similarity  between  the  two  figures,  especially  for 
two  pairs  of  closely  spaced  explosions  (MAST  and  HOSTA,  INLET  and  BACKBEACH). 
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delay-  times. are>snown;in?Figures>20  arid ‘23.  These- figures -include- data5 from -as:  many  as  5 
pairs- of.  closely  spaced:  shots  (denoted  .by  A;B.  and S;D):  for-  the  Yucca1  Flat-region  and  5  pairs 
for  the  P^ihutevMesa  region.  'Nearly.all-closely -spaced  shots  show  differences  in  Pn/Lg  and 
.Pg/Lg  that  agree  well  with ;  the :  theoretically  expected  dependence.  The  results  frojnKNB  are 
therefore, again  ^similar  to  those, from  ELK  (Figures -8 ‘through  :13).  The  similarity  of  results 
from  (ELK  ,and  KNB  (indicates  ;that  the  observed  dependence  of  Pn/Lg  and  Pg/Lg  on  shot 
depth  is  npt.significantlydnfluenced:.by;paih:effects  and  is-mairily  due'to  near-source  scattering 
-of  the  depth-dependent>Rg. 
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COMPARISON  OF  ELK  AND  KNB  RESULTS  FOR  COMMON  SHOTS 


In  order  to  understand  the  influence  of  variations  in  the  propagation  path  on  the  low- 
frequency  spectra  of  regional  phases,  a  comparison  was  made  of  the  spectral  ratios  Pn/Lg  and 
Pg/Lg  derived  from  ELK  and  KNB  records  of  common  Yucca  Flat  shots.  Average  values  of 
Pn/Lg  and  Pg/Lg,  over  the  frequency  range  of  0.5  to  2.0  Hz,  were  plotted  versus  the  estimated 
delay  time  for  22  shots  recorded  at  both  ELK  and  KNB,  and  the  residuals  with  respect  to  the 
least  squares  linear  regression  were  computed.  Residuals  for  various  shots  recorded  at  ELK 
and  KNB  are  plotted  at  their  epicentral  locations  in  Figures  24  and  25  for  Pn/Lg  and  Pg/Lg, 
respectively.  A  comparison  of  the  residuals  at  ELK  and  KNB  for  either  Pn/Lg  or  Pg/Lg  do 
not  show  much  similarity.  In  fact,  many  explosions  show  opposite  trends.  Figure  26  shows 
residuals  at  ELK  plotted  versus  those  at  KNB  for  both  Pn/Lg  and  Pg/Lg.  The  opposite  trend, 
suggested  by  the  negative  mean  slopes,  is  much  better  established  for  Pg/Lg  than  for  Pn/Lg.  It 
is  hard  to  explain  these  large  and  rather  systematic  variations  with  direction  on  the  basis  of 
spall  which  mainly  represents  an  axisymmetric  source.  On  the  other  hand,  theoretical  (e.g. 
Stead  and  Helmberger,  1988)  as  well  as  observational  (e.g.  Gupta  etal.,  1990)  studies  of 
explosion  sources  in  the  Yucca  Flat  region  have  suggested  strong  directional  effects  associated 
with  the  scattering  of  Rg. 
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Figure  24.  Similar  to  Figure  8  for  residuals  of  linear  regression  of  Pn/Lg  (0.5-2.0  Hz)  versus 
estimated  delay  time  derived  from  22  Yucca  Flat  shots  recorded  at  (a)  ELK  and  (b)  KNB. 
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Figure  25.  Similar  to  Figure  8  for  residuals  of  linear  regression  of  Pg/Lg  (0.5-2.0  Hz)  versus 
estimated  delay  time  derived  from  22  Yucca  Flat  shots  recorded  at  (a)  ELK  and  (b)  KNB. 


Figure  26.  Residuals  of  linear  regression  of  (a)  Pn/Lg  and  (b)  Pg/Lg  versus  estimated  delay 
time  derived  from  22  Yucca  Flat  shots  recorded  at  KNB  plotted  against  similar  residuals  from 
ELK.  Note  the  opposite  correlation  (indicated  by  negative  mean  slope)  in  both  cases. 
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SPECTRAL  RATIOS  OF  NORMAL  AND  OVERBURIED  SHOTS 


Figures  27  and  28  show  spectral  ratios  between  two  pairs  of  closely  spaced  normal  and 
overburied  explosions.  BASEBALL  and  BORREGO  (overburied)  have  the  same  shot  depth 
and  are  less  than  500  m  apart  (Figure  1).  The  overburied  shot,  TECH  ADO,  is  only  15  m 
deeper  than  ROUSANNE  and  the  two  are  separated  by  less  than  300  m.  On  the  basis  of  their 
mb  values  (Table  1),  yields  of  the  two  overburied  shots  are  at  least  a  factor  of  about  25 
smaller  than  their  neighboring  explosions.  The  spectral  ratios  are  based  on  the  use  of  S/N  = 
1.5  and  five-point  smoothing.  Results  from  the  three  stations  ELK,  KNB,  and  LAC  are  rather 
similar  and  show  the  low-frequency  Lg  from  the  normal  shots  to  be  significantly  larger  than 
that  from  the  overburied  shots  whereas  such  behavior  is  mostly  not  observed  for  the  other  two 
phases  Pn  and  Pg.  It  is  generally  agreed  that  the  normal  nuclear  shots  are  associated  with  a 
large  non-linear  zone  and  it  is  much  more  realistic  to  regard  them  as  volume  sources  rather 
than  point  sources.  Since  most  of  the  non-linear  zone  is  above  the  shot  point,  the  "effective" 
shot  depth  may  be  considerably  smaller  than  the  emplacement  depth.  On  the  other  hand,  the 
non-linear  zone  associated  with  an  overburied  shot  is  expected  to  be  much  smaller  so  that  its 
effective  shot  depth  may  be  close  to  its  emplacement  depth.  Thus  the  observed  spectral  ratios 
for  the  low-frequency  Lg  in  Figures  27  and  28  can  be  explained  by  the  scattering  of 
explosion-generated  Rg  into  S. 

Theoretical  studies  of  the  generation  of  Lg  by  spall  from  a  normal  explosion  indicate  a 
sharply  peaked  spectrum  (compared  to  the  broadband  explosion  signal)  which  falls-off  rapidly 
for  frequencies  above  2  Hz  and  below  about  0.5  Hz  (Barker  et  al. ,  1990;  McLaughlin,  et  al. , 
1990).  Spall  from  an  overburied  explosion  is  significantly  less  important  than  from  a  normal 
shot  (McLaughlin  et  al. ,  1990).  It  follows  therefore  that  the  amplitude  ratio  of  Lg  for  a 
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Figure  27.  Spectral  ratios  of  closely  spaced  explosions  BASEBALL  and  BORREGO  (over¬ 
buried)  for  phases  Pn,  Pg,  and  Lg  at  stations  ELK,  KNB,  and  LAC.  Note  the  relatively  larger 
low-frequency  Lg  for  BASEBALL  at  all  three  recording  stations. 


ROUSANN  E/TECH  ADO 


normal  to  an  overburied  shot  should  decrease  as  the  frequency  approaches  0  Hz.  Our  obser¬ 
vations  (Figures  27c,  f,  i  and  28c,  f,  i)  clearly  indicate  an  opposite  trend. 
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CONCLUSIONS 


The  observed  variations  of  Pn/Lg  and  Pg/Lg  with  shot  depth  and  delay  time,  especially 
for  closely  spaced  shots  for  which  the  scattering  functions  are  likely  to  be  similar,  suggest  that 
the  low-frequency  Lg  from  Yucca  Flat  and  Pahute  Mesa  (NTS)  explosions  is  mainly  due  to 
the  scattering  of  explosion-generated  Rg  into  S.  The  observed  dependence  of  the  amplitude 
ratios  Pn/Lg  and  Pg/Lg  on  shot  depth  is  mostly  in  agreement  with  the  dependence  on  depth  of 
explosion-generated  Rg.  A  better  quantitative  agreement  with  theory  is  not  expected  because 
of  the  extensive  lateral  variations  in  physical  parameters  that  exist  for  most  NTS  shots.  The 
known  stability  of  Lg  as  an  estimator  of  explosive  yield  may  be  explained  by  scattering  taking 
place  over  a  large  volume  of  the  shallow  crust. 
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SECTION  C 


AN  ARRAY  STUDY  OF  THE  EFFECTS  OF  A  KNOWN  LOCAL  SCATTERER 

ON  REGIONAL  PHASES 

INTRODUCTION 

Until  recently,  very  few  studies  of  teleseismic  PPP  waveforms  had  been  successful  in 
identifying  and  locating  specific  sources  of  scattering  from  observed  seismograms.  One  of  the 
most  convincing  analyses  has  been  the  identification  of  scattered  arrivals  from  the  region  of 
Moffat  Water,  a  narrow  lake  in  Scotland  about  8  miles  long,  which  lies  at  the  bottom  of  a  600 
ft  deep  valley  (Key,  1967).  Broadband  f-k  analysis  of  EKA  recordings  of  both  U.S.  and 
Soviet  nuclear  explosions  further  confirmed  the  secondary  seismic  source  to  be  from  Moffat 
Water  (Gupta  et  al. ,  1990a).  Residual  seismograms  obtained  by  subtracting  the  beamed 
record  from  each  array  channel  (Blandford  etal.,  1976)  significantly  improved  identification 
of  the  sources  of  locally  scattered  arrivals. 

Our  recent  f-k  analyses  of  data  from  the  high-frequency  seismic  array  NORESS  (Nor¬ 
way)  first  suggested  the  presence  of  a  local  scatterer  about  25-30  km  southwest  of  the  array, 
in  the  region  of  Lake  Mjosa  near  Skreikampen  in  Figure  1  (Gupta  et  al.,  1989;  Gupta  et  al., 
1990a,  1990b).  Soviet  and  U.S.  nuclear  explosions  and  an  earthquake  were  used  as  the 
seismic  sources,  providing  significantly  different  source-receiver  azimuthal  directions.  Two 
new  f-k  techniques  based  on  the  use  of  residual  seismograms  and  f-k  power  difference  plots 
improved  identification  of  the  scatterer.  Similar  results  were  obtained  by  polarization  analysis 
of  three-component  NORESS  data.  Our  deterministic  location  of  a  prominent  scatterer  near 
Skreikampen  has  been  duplicated  by  later  studies  (e.g.  Bannister  etal.,  1990;  Hedlin  etal., 


71 


Figure  1.  The  NORESS  array  region  with  the  two  scattering  locations  Skreikampen  and 
Bronkeberget  hatched.  The  array  configuration  is  given  in  the  lower  right  insert  where  encir¬ 
cled  stations  indicate  three-component  sites. 


1991).  We  also  located  another  scatterer  about  10  km  east  of  NORESS,  near  Bronkcberget  in 
Figure  1  (Gupta  et  al. ,  1989)  and  this  too  has  been  confirmed  by  later  research  (Bannister 
ei  al.,  1990;  Hedlin  etal.,  1991).  In  this  study,  we  investigate  the  effects  of  the  well  esta¬ 
blished  prominent  scatterer  associated  with  Lake  Mjosa  on  array  recordings  of  the  regional 
phases  Pn,  Pg,  and  Lg. 

In  order  to  make  best  possible  use  of  array  data,  it  is  important  to  understand  the  charac¬ 
teristics  of  various  arrivals  from  both  teleseismic  and  regional  seismic  sources.  Both  near¬ 
receiver  and  near-source  scattering  play  important  roles  in  the  monitoring  of  underground 
nuclear  tests.  For  example,  the  Lg  phase  from  explosions,  known  to  be  strongly  influenced  by 
scattering  (e.g.,  Gupta  and  Blandford,  1983),  is  often  used  for  the  detection,  source  discrimi¬ 
nation,  and  yield  determination  of  underground  nuclear  explosions. 
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F-K  ANALYSIS  OF  REGIONAL  DATA  FROM  A  MINE  BLAST 


Our  study  of  several  regional  events  recorded  at  NORESS  suggested  that  near-receiver 
scattering  exerts  significant  influence  on  the  regional  phases  Pn,  Pg,  and  Lg  by  generating 
observable  scattered  arrivals  (Gupta  et  al.  1990c).  Results  from  the  mine  blast  of  4  February 
1986  (59.3°  N,  24.4°  E,  mb  =  2.6,  distance  734  km,  backazimuth  97°;  Jurkevics,  1988)  are 
shown  in  Figures  2  and  3.  As  with  teleseismic  data,  the  use  of  residual  seismograms  derived 
by  subtracting  the  beamed  record  from  each  array  channel  provided  an  effective  method  for 
isolating  the  secondary  source  from  the  primary  source.  The  beamed  record  for  each  regional 
phase  was  obtained  by  incorporating  steering  delays  based  on  the  known  azimuthal  direction 
and  phase  velocity  from  the  f-k  analysis  of  raw  data.  Figure  2a  shows  f-k  results  for  the  6.4 
sec  long  window  with  the  largest  amplitude  Lg,  Lgmax.  It  indicates  a  phase  velocity  of  about 
4.6  km/sec  and  a  backazimuth  of  about  96°.  Knowing  the  location  of  the  Skreikampen 
scatterer  and  assuming  a  group  velocity  of  3.5  km/sec  for  Lg,  the  scattered  Lg-to-Rg  arrival 
should  arrive  about  14  sec  after  Lgmax;  the  corresponding  normal  f-k  plot  is  shown  in  Figure 
2b.  F-k  plots  from  residual  seismograms  obtained  by  subtracting  the  beamed  Lg  from  each 
array  channel  provided  the  results  shown  in  Figures  2c  to  2f  in  which  plots  for  both  later  win¬ 
dows  (Figures  2d  and  2f)  eleaily  indicate  the  expected  l.g  to  Rg  arrivals  along  the  ex  pet  ted 
direction  and  with  phase  velocity  appropriate  for  Rg.  A  comparison  of  power  on  the  normal 
and  residual  f-k  plots  suggests  that  amplitude  of  the  Lg-to-Rg  scattered  phase  is  a  significant 
fraction  of  the  on-azimuth  Lg.  For  example,  in  the  Lginax  +  14  sec  window  (Figure  2d),  the 
prominent  scattered  arrival  has  power  only  about  3  db  smaller  than  that  of  the  on-azimuth  Lg 
(Figure  2b). 
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(a)  N.  Lg((inx,  1.25  ±  0.75  Hz 
Az:  96.4  Vol:  4.55  Pwc:  5.9 


(b)  N.  lgni;ix  c  14.0  see.  1.25  +  0.75  Hz 

Az:  130.V  Vcl:  4.05  Pwf:  O.G 


FI.  Uw.  1.25  +  0.75  Hz 

Az:  106.0  Vcl:  13.21  Pwr: -4.2 


(d)  R,  Lgn)ax  +  14.0  sec.  1.25  ±  0.75  Hz 

Az:  249.4  Vd:  2.3C  Pwr: -3.7 


R.  Lgmax,  0.90  ±  0.60  Hz 

Az:  101.1  Vcl:  13.40  Pwr:  -4.7 


(I)  R.  Lg,^  +  14.0  sec,  0.90  ±  0.60  Hz 

A2: 252.7  Vel:  2.42  Pwr:  -4.2 


D 


Figure  2.  Frequency-slowness  spectral  estimates  of  the  mine  blast  of  4  February  1986  with  a 
maximum  slowness  of  0.5  sec/km  and  use  of  Lg  windows.  The  numbers  on  top  of  each  plot 
indicate  backazimuth  (deg),  phase  velocity  (km/sec),  and  signal  power  (dB),  respectively  of 
the  highest-amplitude  arrival  which  is  shown  connected  to  the  plot  origin.  The  signal  win¬ 
dows  and  frequency  passbands  are  indicated  on  each  plot.  Plots  (a)  and  (b)  are  based  on  the 
use  of  normal  (N)  seismograms  whereas  (c)  to  (f)  are  derived  from  the  residual  (R)  seismo¬ 
grams.  Both  (d)  and  (0  indicate  prominent  scattered  Rg  arrivals  from  the  direction  of 
Skreikampen. 


(a) 


0  5  10  If)  20 


TIME  (SRC)  FROM  Lg 

A  NORMAL,  O  RESIDUAL,  1 .25  ±  0.75  Hz 

Figure  3.  Plots  of  (a)  phase  velocity  and  (b)  backazimuth  of  the  dominant  arrivals  in  f-k  plots 
of  25  sec  of  Lg  derived  from  normal  and  residual  seismograms  of  the  mine  blast  of  4  Febru¬ 
ary  1986.  The  size  of  the  symbol  is  proportional  to  the  power  in  each  dominant  arrival. 
Backazimuths  to  the  explosion  and  Slaeikampen  are  denoted  by  EXPL  and  SK,  respectively. 
Arrows  marked  A  and  B  represent  results  from  windows  corresponding  to  Lgmax  and  the 
expected  Lgmax-to-Rg  scattered  phase  from  Skreikampen,  respectively.  Results  from  the  resi¬ 
dual  seismograms  indicate  prominent  Rg  arrivals  from  the  direction  of  Skreikampen. 
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F-k  plots  from  both  normal  and  residual  seismograms  were  obtained  for  25  consecutive 
windows,  starting  5  sec  before  Lgmax  and  with  incremental  shift  of  1  0  sec.  As  before,  each 
window  was  6.4  sec  long.  Results  using  only  those  f-k  plots  which  indicated  the  phase  velo¬ 
city  of  the  dominant  arrival  to  be  within  the  range  of  2.0  to  5.0  km/sec  are  shown  in  Figure  3. 
The  size  of  symbols  is  proportional  to  the  power  in  each  dominant  arrival  but  the  maximum 
power  for  normal  and  residual  f-k  plots  differ  by  about  25  db.  In  Figure  3a,  the  phase  veloci¬ 
ties  of  dominant  arrivals  on  normal  f-k  plots  show  considerable  variation  whereas  the  phase 
velocities  of  most  prominent  arrivals  on  residual  f-k  plots  are  restricted  to  velocities  close  to 
2.5  km/sec,  suggesting  Rg  .  Similarly,  Figure  3b  indicates  nearly  all  dominant  arrivals  on  the 
normal  f-k  plots  arriving  from  directions  close  to  the  backazimuth  to  the  source  and  most 
dominant  arrivals  on  the  residual  f-k  plots  coming  from  the  direction  of  Skreikampen.  The 
arrows  marked  by  the  letters  A  and  B  denote  prominent  arrivals  in  windows  corresponding  to 
Lgmax  a°d  Lgmax  +  14  sec.  The  second  (B)  arrival  is  associated  with  a  phase  velocity 
expected  for  Rg  and  agrees  remarkably  well  with  the  expected  backazimuth  and  arrival  time 
of  Lg-to-Rg  scattered  from  Skreikampen.  The  large  number  of  prominent  arrivals,  with  phase 
velocity  typical  of  Rg  and  azimuthal  direction  close  to  Skreikampen,  observed  on  residual  f-k 
plots  (Figure  3)  suggests  that  a  significant  fraction  of  the  seismic  energy  directly  reaching 
Skreikampen  is  continuously  scattered  into  Rg  recorded  at  NORESS. 
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F-K  ANALYSIS  QF  REGIONAL  DATA  FROM  TWO  LOCAL  EARTHQUAKES 


Results  of  similar  analysis  of  NORESS  data  from  the  earthquake  of  1  April  1986  (56.42° 
N,  12.10°  E,  mj,~  3.6,  distance  482  km,  backazimuth  176°;  Jurkevics,  1988)  are  shown  in 
Figures  4  and  5.  F-k  results  for  the  6.4  sec  long  window  with  the  largest  amplitude  Lg, 
Lgnm,  shown  in  Figure  4a,  indicate  a  phase  velocity  of  about  3.7  km/sec  and  a  backazimuth 
of  about  169°.  Based  on  known  location  of  the  Skreikampen  scatterer,  the  scattered  Lg-to-Rg 
arrival  should  arrive  about  6  sec  after  Lgmax;  the  corresponding  normal  f-k  plot  is  shown  in 
Figure  4b.  F-k  plots  from  residual  seismograms  provided  results  shown  in  Figures  4c  to  4f  in 
which  plots  for  both  later  windows  (Figures  4d  and  4f)  clearly  indicate  the  expected  Lg-to-Rg 
arrivals  along  the  expected  direction  and  with  phase  velocity  appropriate  for  Rg.  In  addition  to 
Lg-to-Rg  arrivals  from  the  direction  of  Skreikampen,  Figures  2d  and  2f  also  suggest  prom¬ 
inent  arrivals  from  a  direction  close  to  that  of  Bronkeberget  and  with  phase  velocity  expected 
of  Rg.  These  secondary  arrivals  are  therefore  probably  due  to  the  scattering  of  Lg-to-Rg  in 
the  vicinity  of  Bronkeberget.  A  comparison  of  the  power  on  the  normal  (Figure  4b)  and  resi¬ 
dual  (Figure  4d)  f-k  plots  again  suggests  that  the  amplitude  of  the  scattered  Lg  is  not 
insignificant. 

Similar  to  the  analysis  of  data  from  the  mine  blast,  f-k  plots  from  both  normal  and  resi¬ 
dual  seismograms  were  obtained  for  25  sec  of  Lg,  starting  5  sec  before  Lgmax.  Results  based 
on  f-k  plots  with  the  phase  velocity  of  the  dominant  arrival  limited  to  2. 0-5.0  km/sec  are 
shown  in  Figure  5.  The  size  of  symbols  is  again  proportional  to  the  power  in  each  dominant 
arrival  but  the  maximum  powers  for  normal  and  residual  f-k  plots  differ  by  about  25  db.  All 
dominant  arrivals  on  the  normal  f-k  plots  are  observed  to  be  from  directions  close  to  the 
backazimuth  to  the  source  whereas  most  dominant  arrivals  on  the  residual  f-k  plots  appear  to 
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AZIMUTH  (DEG)  PHASE  VELOCITY  (KM/SEC) 


(a) 


:  -5  0  5  10  15  20 

TIME  FROM  Lgmax  (SEC) 

A  NORMAL,  O  RESIDUAL,  1.25  ±0.75  Hz 

Figure  5.  Similar  to  Figure  3  for  the  regional  earthquake  of  1  April  1986.  Backazimuths  to 
the  earthquake,  Skreikampen,  and  Bronkeberget  are  denoted  by  EQKJE,  SK,  and  BR,  respec¬ 
tively.  Arrows  marked  A,  B,  and  C  denote  windows  corresponding  to  Lgnm  and  the  expected 
Lgmax-to-Rg  scattered  phases  from  Skreikampen  and  Bronkeberget,  respectively.  Results  from 
the  residual  seismograms  indicate  Rg  arrivals  from  azimuthal  directions  close  to  those  of 
Skreikampen  and  Bronkeberget. 
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be  from  directions  close  to  those  of  the  two  scattering  regions  associated  with  Skreikampen 
and  Bronkeberget.  The  arrows  denoted  by  the  letters  A,  B,  and  C  denote  windows 
corresponding  to  Lgmax,  Lgmax  +  6  sec,  and  Lgmax  +  4  sec,  respectively.  The  last  two  arrival 
times  are  associated  witfT  phase  velocities  expected  for  Rg  and  agree  remarkably  well  with  the 
expected  backazimuth  and  arrival  times  of  Lg-to-Rg  scattered  from  Skreikampen  and 
Bronkeberget,  respectively. 

F-k  plots  were  also  obtained  for  6.4  sec  long  windows  containing  the  Pn,  Pgmax  and  later 
arrivals.  Results  from  normal  seismograms  containing  the  Pn  and  Pgmax  windows  are  shown 
in  Figures  6a  and  6d,  respectively.  The  residual  f-k  plots  for  the  Pn  +  6  sec  and  Pgmax  +  6 
sec  windows  indicated  dominant  arrivals  from  the  south-west  direction  (Figures  6c  and  6f) 
although  the  corresponding  normal  f-k  plots  (Figures  6b  and  6e)  suggested  dominant  arrivals 
directly  from  the  source  region.  Considering  the  phase  velocities  and  the  various  distances 
involved,  it  seems  that  the  primary  arrivals  observed  on  the  two  residual  plots  in  Figures  6c 
and  6f  are  probably  be  due  to  the  scattering  near  Skreikampen  of  Pn-to-Rg  and  Pgmax-to-Rg, 
respectively. 

Similar  analyses  were  carried  on  the  NORESS  data  for  another  earthquake  that  occurred 
along  a  considerably  different  azimuthal  direction.  Results  from  the  earthquake  of  5  February 
1986  (62.81°  N,  4.86°  E,  mb  ~  4.7,  distance  422  km,  backazimuth  306°;  Jurkevics,  1988), 
based  on  the  use  of  Pn  and  Pg  windows,  are  shown  in  Figures  7  and  8,  respectively.  Lg  from 
this  relatively  larger  event  was  clipped  so  that  scattering  of  Lg-to-Rg  could  not  be  investi¬ 
gated.  Pn  from  this  event  was  rather  weak  but  Pg  was  strong  and  impulsive.  In  Figure  7, 
only  the  f-k  plot  with  the  lowest  frequency  passband  (Figure  70  shows  prominent  Rg  arrival 
from  the  direction  of  Skreikampen,  probably  because  of  the  low  amplitude  of  Pn.  However, 
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(a)  N,  Pn  +  1.0  sec,  0.90  ±  0.60  Hz  (b)  N.  Pn  +  6.0  sec,  0.90  ±  0.60  Hz  (C)  H,  nn  +  b.u  sec,  u.yu  ±  u.c 

Az:  194.0  Vet:  6.76  Pwr:  0.1  Az:  161.4  Vel:9.71  Pwr:  2.8  Az:  231.0  Vel:  3.23  Pwr: -4.1 
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Figure  6.  Similar  to  Figure  2  but  based  on  use  of  Pn  and  Pg  windows  on  records  of  the 
regional  earthquake  of  1  April  1986.  Plots  (a),  (b),  (d),  and  (e)  are  obtained  by  using  normal 
(N)  seismograms.  Plots  (c)  and  (f),  derived  from  residual  (R)  seismograms,  indicate  prom¬ 
inent  scattered  arrivals  from  the  direction  of  Skreikampen,  probably  due  to  the  scattering  of 
Pn-to-Rg  and  Pgmax-to-Rg,  respectively. 


(a)  N.  Pn  .  1 .0  see,  1 .25  ±  0.75  Hz  (b)  N,  Pn  +  9  0  see.  1 .25  ±  0.75  Hz 
Az:  304.4  Vel:  0.20  Pwr:  26.4  Az: 307.4  Vcl:G.99  Pwr:  31.9 


(c)  R.  Pn  +  1.0  sec,  1.25  ±  0.75  Hz  (d)  R,  Pn  +  9.0  sec,  1.25  ±  0.75  Hz 
Az:  287.6  Vel:  2.79  Pwr  10.1  Az:  111.3  Vel:  2.76  Pwr:  19.2 


(e)  R.  Pn  4-  1 .0  sec,  0.90  ±  0.60  Hz  (()  R,  Pn  +  9.0  sec,  0.90  ±  0.60  Hz 
Az:  282.5  Vel:  2.30  Pwr  6.0  Az:  257.3  Vel:  3.24  Pwr:  13.9 


Figure  7.  Similar  to  Figure  2  but  based  on  use  of  Pn  windows  on  records  of  the  regional 
earthquake  of  5  February  1986.  Only  the  lowest-frequency  passband  (0  indicates  prominent 
scattered  arrival  from  the  direction  of  Skreikampen,  probably  due  to  the  scattering  of  Pn-to- 
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(e)  R,  Pg**.  0.90  ±  0.60  Hz 

Az:  57.8  Vcl:  17 .35  Pwr:3C.2 


(f)  R.  Pg,,^  +  8.0  sec,  0.90  i  0.60  Hz 
Az:251 .6  Vcl:  2.83  Pwr.36 2 
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Figure  8.  Similar  to  Figure  2  but  based  on  use  of  Pg  windows  on  recotds  of  the  regional 
earthquake  of  5  February  1986.  Both  (d)  and  (0  indicate  prominent  scatteted  Rg  arrivals  from 
the  direction  of  Skreikampen,  probably  due  to  the  scattering  of  Pgmax-to-Rg. 


plots  in  both  Figures  8d  and  8f,  based  on  the  use  of  Pg  windows,  show  clear  evidence  of  Rg 
arrivals  from  the  direction  of  Skreikampen. 

In  addition  to  these  two  events,  several  other  earthquake  and  explosion  sources  along 
various  azimuthal  directions  to  NORESS  also  provided  evidence  for  scattered  arrivals  Pn-to- 
Rg,  Pg-to-Rg,  and  Lg-to-Rg,  originating  mainly  from  the  two  scattering  source  regions  near 
Skreikampen  and  Bronkeberget.  These  scattered  arrivals  arc  observed  with  remarkable  clarity 
and  closeness  to  the  expected  arrival  times  and  are  rich  in  low  frequencies,  as  expected  for 
Rg.  Their  presence  in  various  regional  phases  indicates  that  near-receiver  scattering  plays  an 
important  role  in  the  composition  of  regional  phases. 


85 


POSSIBLE  effect  of  scattering  on  estimation  of  source  azimuth 

Sutcau-Henson  (1990)  examined  azimuth  and  slowness  estimates  by  polarization  analysis 
of  the  first  several  seconds  of  P  arrivals  on  three-component  recordings  of  regional  events 
recorded  at  NORESS.  The  epiccntral  locations  of  these  events  were  known  from  the  Intelli¬ 
gent  Monitoring  System  Bulletins.  Results  for  azimuth  differences  (defined  as  azimuth  from 
polarzation  minus  azimuth  from  the  IMS  Bulletins)  versus  the  IMS  azimuth  for  the  four 
three-component  sensors,  NRAO,  NRC2,  NRC4,  and  NRC7  are  shown  in  Figure  9  (after  Fig¬ 
ure  9b  of  Suteau-Henson,  1990).  Although  NRC7,  which  has  the  largest  number  of  observa¬ 
tions,  shows  considerably  more  deviations  than  the  other  three  sensors,  the  generally  large 
deviation  at  azimuths  of  about  90°  and  230°  -  300°  for  all  sensors  may,  at  least  in  part,  be 
due  to  the  scattering  at  Bronkebergct  and  Skreikampen,  respectively,  which  lie  along  nearly 
the  same  azimuthal  directions  (see  Figure  1). 

Assuming  isotropic  scattering,  the  contribution  of  secondary  arrivals  (such  as  Tn-to-Rg) 
due  to  scattering  at  Bronkeberget  is  expected  to  be  the  earliest  (about  2  sec  after  the  first  Pn) 
for  seismic  sources  lying  east  of  NORESS.  This  may  contaminate  the  Pn  phase  particle 
motion,  producing  larger  errors  in  azimuth  estimates  of  events  east  of  NORESS.  A  similar 
reason  for  larger  azimuthal  errors  may  apply  to  events  lying  along  the  azimuthal  direction  of 
Skreikampen.  A  detailed  study  may  lead  to  azimuth-dependent  corrections  for  the  improve¬ 
ment  of  azimuth  estimates  from  three-component  data  recorded  at  NORESS. 

Odegaard  et  al.  (1990)  investigated  the  effects  of  surface  topography  within  a  few  km  of 
NORESS  and  ARCESS  on  slowness  and  azimuth  anomalies  for  signals  recorded  at  the  two 
arrays.  Surface  topography  was  found  to  explain  about  half  of  the  observed  anomalies  in  the 
three-component  solutions. 
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8  REGIONAL  P  AT  NRAO  25  REGIONAL  P  AT  NRC2 


Figure  9.  Azimuth  differences  versus  theoretical  azimuth  for  regional  P  arrivals  at  the  four 
NORESS  three-component  sensors.  Note  the  generally  larger  differences  along  azimuths 
corresponding  to  Bronkeberget  and  Skreikampen. 


DISCUSSION. 


Use  of  data  from  arrays  such  as  NORESS  is  essential  for  monitoring  decoupled  and 
low-yield  coupled  nuclear  explosions  because  of  the  significant  improvement  in  S/N  provided 
by  multi-channel  data.  Array  measurements  of  Lg  have  also  been  known  to  provide  remark¬ 
ably  precise  measurements  of  yields  of  USSR  nuclear  explosions.  It  is  therefore  important  to 
investigate  methods  which  can  further  improve  measurements  of  regional  phases  at  arrays. 
Our  recently,  developed  methods  including  the  use  of  residual  seismograms  have  already 
demonstrated  their  effectiveness  in  determining  locations  of  prominent  scatterers  near 
NORESS.  A  knowledge  of  the  contribution  of  near-receiver  scattering  associated  with  various 
seismic  sources  (underground  nuclear  explosions  and  earthquakes)  makes  it.  possible  to  correct 
the  observed  Lg  for  the  effects  of  local  (near-receiver)  scattering.  The  improved  estimates  of 
Lg  should  lead  to  improvement  in  Lg-based  yield  estimates  and  source  discrimination. 
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CONCLUSIONS 


F-k  analyses  of  NORESS  array  recordings  of  teleseismic  events  have  indicated  a  prom¬ 
inent  secondary  source  about  25-30  km  southwest  of  the  array,  in  the  region  of  Lake  Mjosa, 
with  large  topographical  relief.  Both  earthquake  and  explosion  sources  at  regional  distances 
and  along  various  azimuthal  directions  to  NORESS  provide  evidence  for  easily  observed  scat¬ 
tered  phases  Pn-to-Rg,  Pg-to-Rg,  and  Lg-to-Rg,  originating  from  the  same  region  of  Lake 
Mjosa.  This  implies  that  scattering  due  to  known  scatterers  near  NORESS  exerts  significant 
influence  on  the  regional  phases  Pn,  Pg,  and  Lg  and  a  detailed  investigation  should  be  useful 
in  improving  location  and  yield  determination  of  regional  events.  Similar  studies  of  the 
influence  of  discrete  local  scatterers  on  regional  phases  should  be  carried  out  for  other  arrays 
such  as  ARCESS  and  GERESS.  The  results  are  useful  not  only  in  the  understanding  of  vari¬ 
ous  arrivals  on  regional  seismograms  but  also  in  improving  yield  estimates  and  source  loca¬ 
tion  and  discrimination. 
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